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[ij  A  combination  of  recent  intensive  observations  and  simulations  with  two  numerical 
models  is  used  to  revisit  the  issue  of  the  northern  Adriatic  response  to  strong  bora 
episodes.  New  observed  and  simulated  data  reinforce  the  view  that  an  episode  of  strong 
bora  wind  provokes  a  double-gyre  (cyclonic,  Trieste,  and  anticyclonic,  Rovinj)  response 
north  of  the  Po  Delta  -  Pula  line.  During  an  intense  bora  episode,  both  measured  and 
modeled  statistics  picture  a  downwind,  highly  polarized,  and  almost  depth-independent 
flow  within  the  Trieste  gyre  NW  ami.  Its  NE  arm  maintains  a  sharp  polarization  and  strong 
depth  dependence  while  exhibiting  lower  speeds,  with  models  in  good  accord  with 
observations.  The  current  statistics  for  Rovinj  gyre  provide  lower  maximum  and  average 
speed  values  and  less  polarized  but  still  rather  depth-independent  flow,  while  exhibiting 
clockwise  rotation.  The  north  arm  of  the  Senj  gyre  (positioned  south  of  the  Po  Delta  Pula 
line)  enjoys  more  lateral  freedom,  and  exhibits  less  rectilinear  (low.  Our  reviewr 
reinforces  the  notion  that  modeling  studies  based  on  ECMWF  wind  forcing  fail  to 
properly  lake  into  account  the  orographic  control  of  the  Dinarie  Alps,  and  to  produce 
correct  bora-induced  gyral  pattern.  The  CO  AMPS  model  successfully  simulated  the 
onset,  duration,  and  decay  of  the  w  ind  peaks,  but  exhibited  a  tendency  to  overpredict  the 
strength  of  the  bora  wind.  Our  simulations  have  identified  the  shallow  NW  coastal  strip 
as  an  important  source  of  colder  water  observed  in  a  sequence  of  remotely  sensed  SST 
fields  derived  from  AVHRR  data. 
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1.  Introduction 

[:]  The  non  hern  Adriatic  is  the  northernmost  part  of  the 
Mediterranean  basin  with  the  Apennines,  Alps,  and  Dinarie 
Alps  providing  its  complex  orographic  surroundings.  One 
of  the  belter  known  winds  there  is  the  bora,  a  cold,  dry,  and 
gusty  downslope  flow  of  air  of  prevailing  northeasterly 
direction.  The  bora  wind  is  an  outcome  of  cold  air  buildup 
over  continental  Europe,  lower  pressure  over  the  Adriatic 
Sea.  and  the  consequent  pressure  gradient  across  the  Dinarie 
Alps,  Although  the  bora  is  closely  related  to  mesoscale 
features  and  synoptic-scale  flow  interactions  with  the 
Alpine  massif  [hancan-Pitek  and  Tutis ,  1996],  gaps  and 
passes  in  the  Dinarie  Alps  (Figure  \)  exert  an  important 
orographic  control,  imposing,  among  other  effects,  lateral 
shearing  in  the  wind  stress  field.  In  a  series  of  numerical 
experiments,  Lazic  and  Tosic  [1998]  showed  that  a  moun¬ 
tain  barrier  about  1000  m  high  is  needed  for  the  occurrence 
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of  bora-type  downslope  wind.  Too  low  a  barrier  prevents 
accumulation  of  cold  air  and  too  high  a  barrier  produces  a 
blocking  effect.  In  one  of  the  first  efforts  to  model  the  ocean 
response  to  lateral  variability  of  the  northern- Adriatic  bora 
wind,  Stravisi  [1977]  used  a  two-dimensional,  storm -surge 
model,  which  assumed  a  sine-squared  function  for  the 
lateral  wind-stress  decay,  Kuzmic  and  Orlic  [1987]  (hereaf¬ 
ter  KG)  used  the  ten-year  statistics  of  Yashina  [1972]  and 
Penza r  [1977]  and  an  interpolation  scheme  to  synthesize  a 
climatological,  laterally  variable,  wind  forcing  field.  When 
applied  to  force  a  three-dimensional,  modal,  numerical 
model  into  a  frictional ly  controlled,  quasi -steady  state,  it 
enforced  a  double-gyre  response  on  the  northern- Adriatic 
circulation  northwest  of  the  Pula-Pesaro  line  (Figure  1 ).  A 
large  cyclonic  gyre  was  obtained  covering  most  of  the 
northern  Adriatic,  with  an  anticyclonic  gyre  apparently 
formed  southeast  of  it.  which  was  dearly  contaminated  by 
the  proximity  of  the  model  open  boundary.  When  the  model 
was  extended  to  cover  the  whole  Adriatic  and  forced  with  a 
similarly  derived  wind  field  [Kuzmic  et  at,  1988],  the 
anticyclonic  gyre  emerged  more  dearly  and  reliably.  This 
work  was  later  expanded  to  include  sirocco  winds  and  more 
complete  interpretation  [Orlic  et  ai.  1994],  The  upwind, 
cross-basin  flow  (part  of  the  cyclonic  gyre)  received  an 
empirical  confirmation  in  a  comparison  with  CZCS  data 
[Kuzmic,  1991]. 
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Figure  1,  The  northern  Adriatic  bathymetry  (contoured,  meters)  and  the  orography  of  the  surrounding 
mountains  (shaded,  meters).  The  three-character  names  mark  the  locations  of  the  ADCP  stations.  Also 
marked  are  the  open  boundary1  of  the  KO  model  (dashed  line)  and  the  names  of  the  locations  mentioned 
in  the  text* 


[}]  In  this  paper,  we  revisit  the  simplistic  KO  simulations 
and  reconsider  them  in  the  light  of  recent  extensive  field 
evidence  and  more  powerful  model  results.  The  purpose  of 
the  revisit  is  threefold:  (a)  to  review  the  modeling  progress 
made  since  the  publication  of  the  KO  paper;  (b)  to  address 
the  validity  of  the  double-gyre  response  with  state-of-the-art 
modeling  tools;  and  (c)  to  verify  the  new  modeling  results 
with  the  currently  available,  much  larger  set  of  field  data. 
The  reanalysis  is  based  on  a  sequence  of  more  sophisticated 
atmospheric  and  oceanographic  model  simulations  and  is 
focused  on  the  January-  February  2003  bora  events*  A  high- 
resolution,  three-dimensional,  finite-element  model  is  used 
to  generate  a  suite  of  simpler,  exploratory  simulations  with  a 
view  to  provide  a  more  insightful  interpretation  of  an  all- 
included*  long-term  simulation  provided  by  another  model. 
In  contrast  to  our  reference  work  (KO)  and  some  other 
studies  dealing  with  (simplified)  generic  processes,  seasonal 
variations,  or  dominant  circulation  patterns,  we  focus  on 
verifiable  predictions  of  several  real  bora  episodes.  The 
novel  points  of  the  revisit  are  the  use  of  two,  sophisticated, 
state-of-the-art  models  of  different  development,  which 
have  (a)  more  realistic  physics,  (b)  higher-resolution  com¬ 
putational  grids,  (c)  hourly  atmospheric  forcing,  which 
allows  detailed  forcing  at  the  air- sea  interface,  and 
(d)  verification  based  on  observations  from  two  recent  field 
programs.  The  improved  verification  relies  on  extensive 


data  sets  collected  at  9  (out  of  I  K)  locations  over  a  period  of 
5  9  months  within  the  framework  of  the  Joint  Research 
Project  (JRP)  (the  U.S.  Naval  Research  Laboratory  (NRL) 
and  the  NATO  Undersea  Research  Center  (NURC))  and  the 
Ruder  Boskovic  Institute  (RBI)  West  ISiria  Experiment 
(WISE).  The  results  of  the  field  measurements  in  the  JRP 
and  WISE  framework,  together  with  other  sources,  provide 
an  unprecedented  empirical  reference  lbr  validation  of  the 
studied  response. 

[4 1  Understanding  of  the  northern  Adriatic  circulation, 
including  observations  and  numerical  simulations,  at  the 
end  of  past  century  has  been  synthesized  in  a  recent  mono¬ 
graph  [Cushman-Roisin  et  aL,  2001],  After  a  relatively  quiet 
period,  the  bora- wind-driven  Adriatic  Sea  circulation  recently 
received  attention  in  several  papers.  Beg  Pakiar  c?  ai  [2001  ] 
performed  a  numerical  model  study  of  the  Adriatic  shelf- 
water  response  to  the  Po  River  discharge  and  a  bora  event. 
The  study  employed  an  MM5  prediction  of  the  wind  stress 
and  surface  heat  flux  to  force  a  version  of  the  Princeton  Ocean 
Model  (POM)  set  up  for  the  Adriatic  at  9- km  resolution.  The 
authors  experimented  with  a  suite  of  forcing  combinations 
that  pointed  out  the  importance  of  proper  treatment  of  air 
stability  and  the  wind- wav  e  field.  Regardless  of  the  treatment 
of  these  parameters,  they  obtained  a  large  cyclonic  gyre  over 
most  of  the  noil  hem  Adriatic,  another  large  cyclonic  gyre 
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south  of' it  (bellow  the  Pula-Pesaro  line.  Figure  1),  and  an 
indication  of  anti-cyclonic  motion  in  between. 

[5]  In  a  four-experiment  numerical  study,  Rachev  ami 
Pur  ini  [200 1  ]  focused  on  the  bora  wind  effects  on  the  Adriatic 
Sea  circulation.  In  these  experiments,  they  forced  an  Adriatic 
implementation  of  the  DieCAST  ocean  model  at  a  spatial 
resolution  of  5.3  km  with  an  idealized  wind  stress  Held  pattern 
based  on  that  used  by  Bergamasco  and  Gacic  [  1  996],  The 
authors  were  primarily  concerned  with  questions  of  the  bora 
contribution  to  the  residual  Adriatic  circulation  and  the 
exchange  at  the  Strait  of  Otranto,  Their  model  grid  resolution 
was  sufficient  to  resolve  details  of  the  bora-induced  circula¬ 
tion  in  the  northern  Adriatic;  however,  the  schematized  nature 
of  their  wind  stress  field  allowed  generation  of  only  a  partially 
realistic  northern  Adriatic  gyral  pattern. 

[fy]  Pullen  et  ai  [2003]  performed  a  realistic,  125-day 
long  simulation  (28  January  4  June  2001)  of  the  Adriatic 
Sea  circulation  with  the  Navy  Coastal  Ocean  Model 
(NCOM)  coupled  to  the  Coupled  Ocean/Atmosphere  Meso- 
scale  Prediction  System  (COAMPS  )  with  a  view  to  identify 
model-derived  patterns  of  circulation  in  both  the  sea  and 
atmosphere.  In  addition  to  assessment  of  the  predicted 
atmospheric  and  oceanic  velocity  fields,  the  authors  set  out 
to  quantify  the  impact  of  the  atmospheric  model  resolution 
on  the  predictive  skill  of  the  ocean  model.  The  double-gyre 
circulation  pattern  was  found  to  be  a  generic  Adriatic 
response  to  the  wind-stress  vorticity  field.  The  high- 
resolution,  4 -km  atmospheric  forcing  produced  more 
realistic  looking  bora  features,  but  low -resolution,  36-km 
model  exhibited  superior  agreement  with  observations  at 
two  A  DC  P- instrumented  sites  for  several  major  statistics, 
(7]  Loglisci  et  ai  [2004]  applied  a  two-way  coupled 
modelling  system  (RAMS-  DieCAST)  to  study  a  severe 
bora  episode  over  the  northern  Adriatic  Sea  that  occurred  in 
January  1995.  Both  models  were  integrated  at  spatial 
resolution  of  7  kin.  The  6-hour  EC  WMF  analyses  of  wind, 
relative  humidity,  and  temperature  were  used  to  provide  the 
initial  and  lateral  boundary  conditions  for  the  atmospheric 
model.  A  elimatologically  initialized,  3 -month-long,  stand¬ 
alone  run  of  the  ocean  model  provided  initial  conditions  for 
the  coupled  system.  Their  results  stress  the  importance  of 
two-way  coupling,  even  for  small  basins  like  the  northern 
Adriatic,  for  quantitatively  correct  simulation  of  thennal 
changes.  However,  the  ECMWF  wind  field  was  too  coarse 
to  depict  the  intricacies  of  the  northern  Adriatic  gyral 
response.  They  obtained  a  very  large  cyclonic  gyre  that 
covered  the  greater  northern  Adriatic.  The  applied  wind 
clearly  missed  the  proper  spatial  structure.  As  a  recent  study 
of  the  Senj  bora  has  shown  [Spoler-Canic  ami  Kraljevic. 
2005]  even  when  the  basic  wind  pattern  is  captured  cor¬ 
rectly,  still  finer  model  resolution  (of  the  order  of  1  km)  is 
likely  to  be  required  to  resolve  the  important  details. 

[hJ  Wang  [2005]  used  a  variant  of  POM  to  investigate  a 
northern  Adriatic  bora  event  in  January  2001 .  ECMWF  wind 
stress  and  heat  flux  fields  were  again  used.  Numerous  (30) 
Adriatic  freshwater  sources  were  represented  as  point  or  line 
source  functions.  The  model  was  run  for  25  months  ( 1  January 
1999  to  31  January  2001)  on  a  regularly  spaced,  5 -km  grid 
with  21  vertical  sigma  levels.  With  such  a  setup,  it  was 
possible  to  simulate  the  weak,  baroclinic.  pre-bora  circulation 
and  changes  triggered  by  the  mid- January  bora.  Like  Loglisci 
et  ai.  [2004],  the  study  primarily  demonstrated  the  bora 


impact  on  the  northern  Adriatic  heat  content  (average  heat 
loss  of  about  200  W/m2),  In  contrast  to  other  previous  studies, 
but  in  accord  with  Loglisci  et  ai  [2004],  Wang  ill  so  predicted 
a  large  cyclonic  gyre  over  the  extended  northern  Adriatic, 
missing  the  small-scale  wind  curl  effect  responsible  for  the 
double-gyre  response.  As  pointed  out  in  a  recent  paper 
[Signed  et  ai.  2005]  the  Adriatic  wind  fields  derived  from 
ECMWF  forecasts  are  not  only  smoothed,  failing  to  properly 
reproduce  the  spatial  structure  of  the  bora  winds,  but  they  also 
significantly  underestimate  the  wind  speed. 

[0]  The  remainder  of  the  paper  is  laid  out  as  follows.  The 
meteorological  and  oceanographic  setup  is  presented  in 
Section  2  and  the  field  data  in  Section  3,  The  mathematical 
models  used  in  the  study  are  introduced  in  Section  4  and  the 
numerical  simulations  elaborated  in  Section  5.  The  results 
are  discussed  in  Section  6  and  the  conclusions  presented  in 
Section  7. 

2,  Meteorological  and  Oceanographic  Setup 

[to]  The  Adriatic  bora  is  a  downs! ope  wind  occurring 
on  the  lee  side  of  the  D inane  Alps,  Its  onset,  longevity, 
and  severity  are  closely  related  10  larger  mesoscate  flows 
[Ivancan-Picek  and  Tads.  1996]  but  ultimately  controlled 
by  local  orographic  features.  Bora  storms  are  connected  to 
two  typical  weather  patterns  [Yoshino,  1976].  The  cyclonic 
pattern  (cyclonic  bora,  “bora  scura")  is  characterized  by  a 
depression  over  the  south  Adriatic,  pulling  the  colder 
continental  air  over  the  mountain  passes.  The  patient  is 
accompanied  by  cloudy  skies  and  occasional  rain.  The 
ant  icy  clonic  pattern  (anticy  clonic  bora,  “bora  chiara")  is 
marked  by  a  high  pressure  system  over  Central  Europe, 
without  necessarily  a  well-developed  cyclone  over  the 
Adriatic,  pushing  the  colder  continental  air  over  the  moun¬ 
tain  passes.  The  pattern  is  accompanied  by  clear  skies,  and 
dry  air.  In  either  case,  accumulation  of  cold  air  landward  of 
the  Dinaric  Alps  is  essential  as  well  as  development  a  strong 
pressure  gradient,  across  the  mountain  gaps  in  particular. 
The  gaps  exert  orographic  control  over  the  air  flow, 
producing  lateral  vorticity  and  “fingers"  of  higher  wind 
speed  (stress)  leeward  and  offshore  of  the  mountain  range. 
Upstream  accumulation  of  cold  air  is  more  likely  to  take 
place  during  the  colder  part  of  the  year,  so  bora  events 
occur  more  often  during  the  winter. 

[11]  Northern  Adriatic  oceanographic  conditions  are 
strongly  influenced  by  its  physiographic  features:  the  shal¬ 
lowness  of  the  water  column,  the  moderate  bathymetric 
gradient,  and  the  sizable  riverine  inflow  [Franco  and 
Michelato.  1992],  Progressive  heat  loss  during  the  winter 
combined  with  wind-induced  mechanical  mixing  destabilizes 
the  water  column  and  enforces  almost  complete  vertical 
homogeneity.  Lateral  temperature  and  salinity  variations 
persist  throughout  the  winter  period,  primarily  along  the 
western  coast,  maintained  by  the  freshwater  inflow  from  the 
Pound  other  rivers.  This  view  is  reinforced  by  a  recent  review 
[Pau lain  etal. ,  2001  ]  ,  and  recently  compiled  temperature  and 
salinity  climatology  (I.  Janekovie,  internal  report,  2005). 

3.  Data  Collection  and  Reduction 

[d]  The  January -February  2003  period  considered  in 
this  study  is  characterized  by  several  strong  bora  events. 
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Figure  2.  Stick  plot  of  the  10-m,  hourly  wind  registered  at  five  east-coast  stations  during  the  January 
February  2003  period.  The  dashed  bars  mark  the  beginning  and  end  of  the  two  prolonged  bora  episodes 
selected  for  study. 


The  January  events  were  distinguished  by  pronounced 
cvclogenesis  over  the  Tyrrhenian  Sea  and  intense  anticy¬ 
clones  over  northwestern  Europe  imparting  strong  pressure 
gradients  over  the  Dinar ic  Alps  (B,  IvancarvPicek,  personal 


communication)*  The  mid- February  bora  event  was  pro¬ 
voked  by  a  strong  anticyclone  north  and  northeast  of  ihe 
Alps.  The  January  events  fit  the  cyclonic  bora  classification 
by  notable  lack  of  AVHRR  scenes  (cloudiness).  During  the 


4  of  27 


C03S13 


KUZM1C  ET  AL:  NORTI U  RN  ADRIATIC  RESPONSE  TO  BORA  WIND 


C03SI3 


mentioned  period  the  10-m  wind  was  registered  at  several 
meteorological  stations. 

[u]  A  stick  plot  of  hourly  wind  vectors  (oceanographic 
convention)  registered  during  the  analysis  period  at  5 
eastern  Adriatic  meteorological  stations  {Trieste,  Rovinj, 
Pula,  Senj,  Mali  Losinj  -  Figure  l )  is  shown  in  Figure  2, 
The  position  of  the  Senj  station  has  been  suspected  to  cause 
under- rep  resen  tat  ion  of  strong  bora  wind  magnitude,  while 
correctly  registering  the  wind  direction.  In  response  to  that 
concern,  we  have  applied  to  the  Senj  magnitude  recently 
derived  correction  factors,  based  on  a  comparison  of  wind 
magnitude  at  the  regular  and  an  auxiliary  (better  positioned) 
Senj  station  (Z.  Bencetic  Klaic,  personal  communication). 
The  correction  increased  the  wind  magnitudes  by  07,  2,7, 
5.0,  and  6.K  ms  1  in  the  0-5,  6-10,  11-15,  and  >15  ms  1 
magnitude  ranges,  respectively 

[u]  One  readily  observ  es  the  prominence  of  northeasterly 
winds  and  their  intermittent  nature.  Closer  inspection  of  the 
figure  reveals  that  the  strongest  wind  was  registered  at  the 
Senj  station,  whereas  the  Rovinj  station,  situated  between 
the  Trieste  and  Senj  bora  corridors,  registered  bora  events  of 
much  lower  intensity.  One  also  notes  the  simultaneity  of  the 
Trieste  and  Senj  bora  events,  but  readily  observes  that  they 
are  not  strictly  coincident  or  of  the  same  intensity  or 
duration.  If  one  selects  the  Senj  record,  takes  the  ENE 
and  NE  as  bora  directions,  and  sets  5  ms  1  as  a  threshold 
speed,  tw'o  periods  of  unusually  long  bora  duration  emerge: 
5-13  January  and  9  17  February.  If  the  speed  threshold  is 
lifted  to  10  ms  ',  a  somewhat  shorter  January  period  arises 
(7  13),  while  a  gap  of  lower  intensity  bora  appears  in  the 
other  episode  on  15  February,  followed  by  two  more  bora 
events.  The  bora  events  in  the  Trieste  time  series  appear 
more  intermittent  and  of  shorter  duration  compared  to  the 
Senj  data.  In  addition  to  the  stated  magnitude  and  direction, 
we  have  required  the  parameters  to  apply  for  at  least  3  days 
and  be  as  applicable  as  possible  to  the  Trieste  record  as  well. 
With  these  requirements,  we  ended  up  with  two  represen* 
tative  bora  events,  one  for  January  (7-  13)  and  another  for 
February  (9  14)  for  closer  modeling  study  and  analysis. 
Both  the  Senj  and  Trieste  series  show  several  additional 
intense  bora  events  of  shorter  duration  and  of  lowrer  Senj* 
Trieste  simultaneity.  It  is  understood  that  alternative  bora 
event  definitions  are  possible,  but  the  one  taken  appears  to 
serve  well  the  stated  purpose  of  the  present  study.  For 
example,  in  their  paper  on  marine  atmospheric  conditions 
and  bora  over  the  nonhem  Adriatic,  Dorman  et  al.  [2006] 
declare  existence  of  bora  at  their  reference  station  when  the 
wind  is  blowing  from  60°  (true  north)  with  a  speed  greater 
than  2,6  ms  1  longer  than  24  hours. 

[^\  Intensive  international,  multidisciplinary  Held  studies 
undertaken  in  the  northern  Adriatic  during  winter  2002/03 
documented  its  oceanographic  conditions  with  an  unprece¬ 
dented  level  of  observational  intensity  [Lee  et  al. ,  2005a], 
NRL  and  NURC  jointly  deployed  14  upward  looking, 
bottom-mounted.  Acoustic  Doppler  Current  Profilers 
(ADCPs)  [Perkins  et  al ..  2000]  distributed  along  3  transects 
in  the  nonhem  Adriatic  and  one  across  Kvamer  Bay.  An 
additional  upward- looking  ADCP  was  mounted  near  the 
base  of  research  tower  Aqua  Alta  (VR1  in  Figure  1).  The 
circle-marked  stations  in  Figure  1  have  been  used  in  present 
study.  At  the  locations  of  the  JRP  moorings  measurements 
were  made  of  currents  throughout  the  water  column  (by 


Table  I .  J  R  P  and  W  l  S  E  ADC  P  Dep  loy  m  en  l  De  lai  I  s'1 


Station 

Longitude 

Latitude 

Depth 

Cell 

Observation  Period 

VR1 

12.51 

45.31 

15.2 

0.5 

04,09.2002, 

0406.2003. 

VR2 

12.64 

45.28 

21.4 

l 

08,10.2002. 

04,05.2003. 

VR4 

1 3.02 

45.18 

30,6 

0.5 

24.09.2002. 

29  >04.2003. 

VR5 

13.28 

45.12 

31.2 

l 

24,09.2002 

04  05.2003. 

VR6 

13,54 

45.06 

29.6 

1 

24.09.2002. 

04.05.2003. 

IC 

13.59 

44.85 

41.2 

2 

20.12.2002. 

07.03.2003. 

CP5 

13,62 

44.68 

43,9 

2 

20.12.2002. 

07.03.2003. 

CP2 

12,86 

44.46 

38.4 

1 

23.09.2002. 

28.04.2003. 

CP3 

13.13 

44.54 

38.4 

1 

23.09.2002 

04  05.2003. 

“Stations  coordinates  arc  in  decrees;  location  and  ceil  depths  are  in 
meters. 


ADCP),  bottom  temperature  (by  ADCP  and  at  some  sites  by 
wave/tide  gauge),  and  bottom  pressure  (by  ADCP  or  wave/ 
tide  gauge).  To  improve  near-surface  results  by  reducing 
surface  wave  aliasing,  the  JRP  ADCPs  w^ere  set  to  measure 
the  currents  using  bursts  of  pings  every  15  minutes  at  I -Hz 
sampling  frequency,  f  urther  processing  details  are  given  by 
J.  W.  Book  et  al.  (Measurements  of  storm  and  nonstorm 
circulation  in  the  northern  Adriatic:  October  2002  to  April 
2003,  submitted  to  Journal  of  Geophysical  Research.  2006, 
hereinafter  referred  to  as  Book  el  al,  submitted  manuscript, 
2006). 

We  have  to  point  out  here  an  anomalous  result 
obtained  for  the  VR5  station.  Extensive  ensuing  analyses 
have  demonstrated  considerable  directional  discrepancy  at 
VR5,  at  odds  with  results  for  the  neighboring  stations  VR4 
and  VR6.  Although  no  inti icat ion  of  instrumentation  error 
was  found,  two  separate  tidal  studies,  one  using  an  incre¬ 
mental-assimilation-augmented,  very-high  resolution  Adri¬ 
atic  tidal  model  [Janekovic  et  al !.,  2004],  and  the  other  a 
strongly  constrained  assimilation  of  JRP  current  and  pres¬ 
sure  data  (Book  et  at,,  submitted  manuscript,  2006)  have 
suggested  the  need  for  a  directional  clockwise  correction  of 
the  VR5  tidal  ellipses  of  25°  and  28°,  respectively.  Such  a 
correction  is  also  in  accord  with  the  bora  field  and  modeling 
results  of  the  present  study.  We  have  therefore  decided  to 
apply  a  directional  correction  to  the  VR5  data,  opting  for 
more  conservative  25°  value, 

[r ?]  The  RBI  WISE  was  conceived  with  a  view  to 
improve  understanding  of  the  northeastern  Adriatic 
response  to  intense  bora  episodes  and  was  partly  modified 
at  the  final  stage  to  fit  the  JRP  mooring  layout.  The  field 
component  comprised  deployment  of  three  ADCPs  in  trawl- 
resistant  bottom  mounts  (two  are  used  in  the  present  study: 
JC  and  CP5  in  Figure  l).  The  WISE  ADCPs  used  uniform 
1 2-sec  time  gaps  between  pings,  2-m  cells,  and  10-min 
sampling  periods.  A  raw  data  processing  similar  to  the  JRP 
procedure  was  applied  to  the  WISE  data,  including  the  same 
response-method-based  tidal  signal  removal  and  2-hour 
low-pass  data  smoother.  However,  the  surface  echo  inter¬ 
ference  was  handled  in  a  more  restrictive  fashion  by 
excluding  the  contaminated  surface  bins  for  the  whole 
measurement  period.  Mooring  details  for  both  field  experi¬ 
ments  are  summarized  in  Table  l .  Both  the  JRP  and  WISE 
data  were  further  low -pass- filtered  w  ith  an  Sth-order  Butter- 
worth  filter  with  a  cutoff  period  of  24  hours  (tides  were 
already  removed  by  the  response  method)  to  remove 
remaining  shorter- period  variance  that  we  have  not  aspired 
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Ta ble  2 a.  L o w -Pass- F i I ic red  C u rren I  Statistics:  VRla 


Period 

Depth 

W 

W 

iii-ai. 

ffT 

0 

u 

vMt 

V 

V^in 

Jan- 

1.7 

14.1 

36.7 

79  1 

54,8 

12.9 

9.7 

1.6 

3,8 

0.39 

Feb 

6.2 

10.4 

31.9 

594 

50. ! 

9.3 

8,5 

1.0 

2,6 

0*30 

15,2 

7.5 

29.1 

46.8 

47.6 

5.8 

K.l 

-02 

2,1 

0,26 

vbar 

9.4 

30*6 

56.7 

48.7 

8.6 

8.2 

0.4 

2*0 

0.25 

0,87 

Jan 

1.7 

16.3 

31*7 

82.8 

55,8 

15,8 

9.0 

2,3 

1.8 

0.19 

episode 

6,2 

15.7 

31*9 

78.7 

55,4 

15*1 

9.4 

2,6 

1.6 

0*17 

15*2 

13*2 

27.0 

62.8 

46*5 

13.0 

8,1 

0.1 

1.7 

0.21 

vbar 

15. 1 

10*3 

75, H 

54*1 

14.7 

9.0 

2,0 

1.5 

0  17 

0,98 

Feb 

1.7 

20,0 

30.5 

35,7 

52*3 

19.8 

6*0 

2,6 

1,2 

0*20 

episode 

6,2 

18.6 

30.3 

39,6 

59.8 

17,8 

6.5 

4.6 

1.4 

0.21 

15,2 

14.4 

26.7 

50.6 

48.8 

14,1 

7.4 

-0,4 

1.9 

0.25 

vbar 

17.5 

29,5 

43.1 

57.6 

r.  i 

6.7 

3.5 

1.2 

0.18 

0.93 

■'W*  mean  speed;  \Vimv*  maximum  speed:  rrT.  total  variance;  il  direction  of  the  major  axis,  li*  mean  current  along  major  axis;  v  rrM.  standard  deviation 
along  major  axis;  v\  mean  current  along  minor  axis;  y<rmt  standard  deviation  along  minor  axis:  :  =  ■  polarization  of  the  flow;  ■ .  current  depth 

dependence:  vbar  marks  the  column  averaged  values.  Depth  is  in  meters*  speed  is  in  ems  angle  in  degrees  north 


to  reproduce  with  the  models.  Current  statistics  fora  three- 
station  subset  of  the  nine  JRP/WISE  stations  used  in  the 
present  study  are  summarized  in  Tables  2a  2c  for  the  whole 
January  February  period  as  well  as  for  the  selected  bora 
events.  The  calculated  statistics  includes  the  mean  and 
maximum  speed  ( W  and  £Knvax),  total  variance  («r7),  orien¬ 
tation  of  the  principal  axis  of  variability*  mean  currents  and 
their  standard  deviations  along  major  and  minor  principal 
axes  (0t  ti*  y/a\u  v*  and  y/<rm%  respectively),  polarization  (e), 
and  current-depth  dependence  (~.)  calculated  following  the 
work  of  Thompson  and  Pugh  [1986]  as; 

(i) 

where  N  is  the  number  of  depth  cells  used  to  form  the  depth 
mean,  i  is  the  index  of  the  depth  cell  summation,  and  the 
overbar  indicates  the  principal  variances  of  the  depth- 
mean  current  Principal  component  analysis  [Emery  and 
Thompson t  1997]  was  used  to  calculated  the  values  of  rrjV/, 
ann  and  C-).  The  selected  stations  (VRI*  IC,  and  CP5) 
demonstrate  the  dominant  response*  At  VR1  (more  influ¬ 
enced  by  the  Trieste  bora),  one  can  readily  observe  that  the 
time*  and  depth-averaged  speed  is  twice  as  large  during  the 
bora  episodes  than  for  the  whole  period.  For  both  episodes, 
the  axis  of  principal  variance  follows  the  local  coastline  to 
within  3  4°,  but  some  other  parameters  exhibit  inter- 
episodic  differences.  Compared  to  February,  the  January  r 
and  7  values  suggest  more  polarized  and  more  depth- 


independent  motion*  but  smaller  averaged  velocity  along  the 
axis  of  principal  variance*  The  case  is  opposite  at  CP5,  a 
station  more  sensitive  to  the  Senj  bora,  where  the  January 
episode  is  characterized  by  less  rectilinear  and  more  depth- 
dependent  flow*  There  is  also  a  factor  of  2.5  difference  in  the 
inter-episodic  orientation  of  the  principal  axis.  The  IC  station 
is  in  the  wind-quiet  region  between  the  Trieste  and  Sen  j  bora 
corridors  with  clearly  lower  average  speeds  and  less 
polarized  flow.  The  CP2  data  for  both  episodes  (not 
tabulated)  exhibit  even  lower  average  and  maximum  speed 
values*  Tabulated  surface,  middepth,  and  bottom  values 
suggest  that  during  the  winter  of  2003  and  the  bora  episodes 
in  particular,  the  northern  Adriatic  water  column  can  be 
viewed  as  a  single  layer.  We  will  return  to  these  statistics  in 
later  sections  when  analyzing  the  modeling  results, 

[is]  A  scatterplot  of  the  vertically  averaged  currents  for  all 
nine  stations  and  the  two  bora  episodes*  given  in  Figure  3* 
provide  additional  insight.  For  both  episodes,  live  self¬ 
similar  groups  can  be  identified  in  the  figure.  The  VR1 
and  VR2  plots  show  downwind  transports  whereas  stations 
VR4  and  VR5  show  a  return*  upwind  flow*  The  third  group 
(VR6  and  IC)  suggests  anticyclonic  motion  along  the  west 
Istrian  coast*  The  CP  stations  3  and  5  indicate  the  existence 
of  another  cyclonic  cell,  while  CP2  appears  outside  any  of 
the  mentioned  circulation  structures. 

4.  Numerical  Models 

[w]  Three  numerical  models  were  used  in  this  study*  one 
atmospheric*  and  two  oceanic.  The  atmospheric  portion  of 


Table  2b.  Low-Pass-Filtered  Current  Statistics:  ICil 


Period 

Depth 

W 

W™. 

£l 

i) 

u 

v/^M 

V 

v 

Jan- 

7*2 

5.3 

16.7 

115 

-65.6 

—2.2 

5.4 

-0*2 

2.5 

0.46 

Fcb 

19,2 

4.8 

15.6 

10.9 

-63.8 

-2,4 

4,9 

0.4 

2.1 

0.43 

39.2 

4,4 

13.1 

7*7 

-52,6 

-3.0 

3.8 

0.4 

2.0 

0.54 

vbar 

4,6 

14,8 

9*8 

-61.6 

-2.5 

4.5 

0.4 

2.0 

0.43 

0.93 

Jan 

7,2 

9,6 

16.7 

19.5 

-80,8 

-8.8 

4.9 

2.6 

2,1 

0.43 

episode 

19.2 

9*3 

15,6 

17,4 

-85.2 

-7.7 

4.9 

3.6 

2.5 

0.5  1 

39.2 

7,9 

13.1 

10.0 

84.6 

-5.6 

3.8 

4.3 

2,6 

0.68 

vbar 

8.9 

14.9 

14.9 

-87,9 

-7,2 

4.5 

3.8 

2,4 

0.54 

0.97 

Feb 

7,2 

49 

8.2 

4.1 

-88,7 

-2.8 

4.1 

-1.0 

1,8 

0.43 

episode 

19.2 

4,5 

7.4 

3.4 

-  80,8 

-2.4 

3,9 

-0.1 

17 

0,44 

39,2 

4.1 

7J 

4,5 

-75,0 

-3.0 

3,1 

0*6 

1.7 

0.54 

vbar 

4.3 

7.3 

3.8 

79.6 

-2.6 

3*6 

0,1 

1*5 

0.42 

0.94 

3 As  in  Table  2a. 
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Table  2c.  Lo  w  -  Pass- Hke  red  Current  Statistics:  CP5fl 


Period 

Depth 

W 

^  IllJX 

£l 

$ 

u 

V 

iMu 

-j 

Jan- 

7.9 

10.1 

30 

62.0 

11.8 

8.5 

8*8 

“0,6 

3.8 

0,43 

Feb 

21.9 

9.1 

36,6 

56.2 

1  u 

7.7 

8.2 

-0,4 

3.5 

0,43 

39.9 

8.2 

33.0 

45.0 

14.5 

6.7 

7,4 

14 

3,5 

0,47 

vbar 

8.8 

35.2 

53.2 

111 

7.4 

7.9 

-0,8 

3.4 

0,42 

0.96 

Jan 

7.9 

19.4 

38.1 

99.0 

6.1 

17.6 

!  1.5 

-2.9 

4.9 

0.43 

episode 

21.9 

17.8 

36.3 

97.6 

6.6 

16.4 

10.8 

-23 

4.9 

0.45 

39.9 

15.4 

33.1 

86.3 

116 

14,3 

9.8 

-2.0 

4.3 

0.44 

vbar 

17.2 

35,3 

94.5 

7.6 

159 

10,5 

—2.4 

4.5 

0.42 

0.96 

Feb 

7.9 

12.4 

23,5 

38.0 

19.0 

12.1 

6.3 

1.5 

1,9 

030 

episode 

21.9 

11.6 

22,8 

36.0 

19.1 

11.3 

6.2 

1.9 

1.4 

0,23 

39.9 

10,6 

20,5 

30.7 

29.8 

10,2 

5.7 

2.3 

1.4 

0.25 

vbar 

(1  3 

21,7 

34.5 

21  8 

11,0 

6,0 

19 

L2 

0,19 

0,96 

'As  in  Table  2a. 
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Figure  3.  Scatterplot  of  the  low -pass-filtered,  vertically  averaged  currents  measured  at  nine  stations 
during  the  January  and  February  bora  episodes. 
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ami: 


the  COAMPS  [Hodun  I W7]  was  applied  in  a  reanalysis 
mode  in  order  to  generate  atmospheric  forcing  fields  needed 
to  drive  the  oceanic  models.  The  COAMPS  atmospheric 
model  is  a  finite -difference  approximation  to  the  fully 
compressible,  non-hydrostatic  equations.  Physical  parame- 
terizations  are  used  to  represent  surface  fluxes,  boundary 
layer  mixing,  radiation,  and  moist  processes  including 
microphysical  quantities  [Hodun  1997],  The  initial  fields 
for  the  model  are  created  from  multivariate  optimum 
interpolation  analyses  of  upper-air  soundings  and  surface, 
commercial  aircraft,  and  satellite  data  that  are  quality 
controlled  and  blended  with  the  12-h  COAMPS  forecast 
fields.  The  data  assimilation  is  accomplished  through  an 
incremental  update  procedure  that  enables  mesoscale  phe¬ 
nomena  to  be  retained  in  the  analysis  increment  fields.  The 
lateral  boundary  conditions  for  the  outermost  mesh  make 
use  of  Navy  Operational  Global  Analysis  and  Prediction 
System  (NOGAPS)  forecast  fields.  The  domain  configura¬ 
tion  for  these  reanalysis  simulations  contains  three  horizon¬ 
tally  nested  grid  meshes  with  horizontal  grid  resolution  of 
36,  12,  and  4  km,  respectively.  The  4-km  resolution  grid 
mesh  is  centered  over  the  Adriatic  Sea.  The  model  is 
configured  with  30  vertical  levels  on  a  non-uniform  vertical 
grid  consisting  of  an  increment  of  10  m  at  the  lowest  level. 
The  topographic  data  for  the  atmospheric  simulations  are 
based  on  the  National  Imagery  and  Mapping  Agency 
(NIMA)  l -km  resolution  data  set  that  enables  the  prominent 
orographic  features,  and  the  Dinaric  Alps  in  particular,  to  be 
resolved.  Additional  details  of  the  COAMPS  reanalysis  can 
be  found  in  Pullen  et  al.  [2003]. 

[20]  The  first  ocean  model  used  here  is  NCOM  as 
described  by  Martin  [2000],  with  some  improvements  as 
reported  by  Morey  et  al.  [2003]  and  Barron  et  al.  [2005], 
This  model  is  similar  in  its  physics  and  numerics  to  POM 
[Blumherg  and  Mellon  1987],  but  uses  an  implicit  treatment 
of  the  free  surface  and  a  hybrid  vertical  grid  with  sigma 
coordinates  in  the  upper  layers  and  (optionally)  level 
coordinates  below  a  user-specified  depth.  For  the  results 
presented  here,  a  third-order  upwind  scheme  [Holland  et  al., 

1 998]  was  used  for  advection,  vertical  mixing  was  computed 
using  the  Mellor-Yamada  Level  2  scheme  [Me! lor  and 
Yam  a  da,  1974],  which  is  modified  for  use  over  the  entire 
water  column,  and  density  was  computed  using  the  equation 
of  state  of  Me l tor  [1991],  The  model  domain  consists  of  the 
entire  Adriatic  Sea  including  the  Strait  of  Otranto.  The 
horizontal  grid  resolution  is  1  km.  The  vertical  coordinate 
consists  of  32  layers,  with  22  sigma  layers  between  the 
surface  and  a  depth  of  290  m  and  level  coordinates  below 
290  m.  Hence,  the  grid  is  like  a  regular  sigma  coordinate 
grid  in  water  shallower  than  290  m  and  is  similar  to  a  level 
grid  where  the  bottom  depth  is  greater  than  290  m.  The 
static  thickness  of  the  upper  layer  is  2  m  w  here  the  bottom 
depth  exceeds  290  m  and  is  reduced  in  shallower  whaler 
where  the  sigma-coordinate  grid  is  compressed.  The  time 
step  was  200  s, 

[zi]  The  second  oceanic  model  is  a  variant  of  the  finite 
element  model  QUODDY  with  some  improvements  based 
on  the  3D,  non-linear,  shallow- water  equations  [Lynch  et 
al. ,  1996].  The  Level  2,5  turbulence-closure  scheme  of 
Mel  hr  and  Yam  ad  a  [1982]  is  used  with  improvements 
described  by  Galperin  et  al.  [1988].  The  horizontal  diffu¬ 
sion  parameterization  scheme  is  that  of  Smagoritisky  [1963], 


A  Tree-slip  condition  is  imposed  along  the  coast.  The  details 
of  the  model  solution  on  a  finite  element  grid  are  described 
by  Lynch  and  Werner  [1991].  Bottom  stress  is  estimated 
from  the  classical  quadratic  law  as  a  function  of  bottom 
velocity.  The  well-known  sigma-coordinate  system  is  used 
in  the  vertical  with  21  non-uniformly  placed  nodes  whose 
sinusoidal  vertical  spacing  provides  increased  resolution  in 
the  surface  and  bottom  layers.  The  unstructured  model  mesh 
consists  of  28669  elements  and  17284  nodes  with  minimal 
nodal  distances  found  in  the  coastal  areas  of  about  100  m 
and  minimal  element  size  about  9000  m\  Maximal  nodal 
distances  of  5500  m  are  found  in  the  middle  of  the  domain. 
With  the  southern  open  boundary  set  at  43.5CN.  the  model 
domain  covers  part  of  the  central  and  the  entire  northern 
Adriatic.  In  accord  with  the  CFL  criterion,  the  model  time 
step  was  set  to  5  s.  Atmospheric  forcing  was  provided  by 
hourly  COAMPS  re-analysis  H)-m  wind  fields. 

5*  Computer  Simulations 

[22]  A  sequence  of  numerical  experiments  was  conducted 
in  order  to  re-examine  the  question  of  the  double-gyre 
response.  Simplified,  team  re- oriented  numerical  simulations 
w'ere  carried  out  using  the  QUODDY  model  whereas  the 
reference  bora  event  simulations  were  taken  from  the  long, 
“everything  included'’  NCOM  Adriatic  run.  In  accord  with 
previous  experience,  the  episodic,  bora-event  QUODDY 
simulations  were  initialized  from  rest  with  a  spin  up  time  of 
several  days.  A  value  of  0.003  was  used  for  the  bottom  drag 
coefficient,  based  on  previously  performed  numerical 
experiments  and  reported  literature  values.  At  the  surface 
the  drag  coefficient  was  defined  according  to  Large  and 
Pond  [1981].  The  river  inflows  were  not  simulated.  Tidal 
forcing  was  not  included  in  the  QUODDY  simulations, 
although  a  separate  tidal  study  w-as  made  [Janekovic  and 
Kuzmin  2005]. 

[23]  The  NCOM  Adriatic  simulation  was  initialized  on 
1  September  2002  with  fields  from  a  hindeasi  ofNRL's  global 
NCOM  model  [Kara  et  at. ,  2005],  Boundary  conditions  at 
the  open  boundary  in  the  northern  Ionian  Sea  were  provided 
by  daily  values  of  surface  elevation,  velocity,  temperature, 
and  salinity  from  the  global  model.  Tidal  forcing  was 
provided  at  the  open  boundaries  using  tidal  data  from  the 
Oregon  State  University  (OSU)  tidal  databases,  which  are 
derived  from  satellite  altimetry  data  [Egbert  and  Erofeeva, 
2002],  Data  from  the  OSU  regional  Mediterranean  database 
were  used  tor  the  KL  01,  M2,  and  S2  constituents  and  data 
from  the  OSU  global  database  were  used  tor  PL  QL  K2, 
and  N2.  Tidal  potential  forcing  for  these  eight  constituents 
was  used  in  the  interior  of  the  model  domain.  Atmospheric 
forcing  was  provided  by  hourly  COAMPS  re-analysis  fields 
of  surface  air  pressure,  wind  stress,  solar  radiation,  net 
longwave  radiation,  and  precipitation  [Hod it n  1997].  Latent 
and  sensible  heat  fluxes  were  computed  using  standard  bulk 
formulas,  the  COAMPS  10-m  wind  speed  and  2-m  air 
temperature  and  humidity,  and  the  NCOM  model  SSL 
The  stability-dependent  Kondo  [1975]  drag  coefficient 
was  used  for  the  bulk  flux  calculations  with  neutral  values 
of  0.0014  and  0.0011  for  the  latent  and  sensible  fluxes, 
respectively.  River  and  runoff  inflows  for  the  Adriatic  were 
taken  from  the  monthly  climatological  database  of  Raieirh 
[1994].  The  discharges  for  39  Adriatic  rivers  were  specified 
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at  the  closest  coastal  model  grid  point  to  the  river-mouth 
location  and  runoff  inflows  were  distributed  along  the 
appropriate  sections  of  coastline,  flatly  observed  discharge 
values  were  used  only  for  the  Po  River  (obtained  from 
Alessandro  Allodi  at  the  Ufficio  Idrografico  del  Magistrate 
per  il  Po,  Parma,  Italy).  The  model  was  run  to  the  end  of 
March  2003  and  model  fields  were  saved  every  3  h  for 
analysis  and  comparison  with  observations, 

[24]  Numerical  experiments  were  conducted  to  address 
three  related  response  issues:  the  time- variable  bora  wind, 
the  gyral-response  pattern,  and  the  extrusion  of  colder/ 
freshened  waters.  The  KG  paper  dealt  with  these  issues 
under  simplified  and  therefore  restrictive  assumptions.  The 
key  issue  of  this  revisit  is  thus  the  wind-imparted,  double¬ 
gyre  pattern  with  additional  aspects  that  accompany  bora 
events:  the  cross-basin  protrusion  of  freshened  water  and 
pronounced  heat  loss  at  the  air-sea  interface.  The  issues  are 
directly  related  to  the  state  of  the  atmosphere  above  the 
northern  Adriatic  and  its  effect  imparted  on  the  NCOM 
and  QUGDDY  model  solutions  via  the  COAM  PS-derived 
forcing, 

[is]  To  address  the  first  issue,  the  COAMPS  10-m  wind 
field  time  series  for  the  two  selected  bora  episodes  were 
compared  to  wind  measurements  at  five  coastal  stations 
(Trieste,  Rovinj,  Pula,  Mali  Losinj,  and  Senj).  Both  models 
were  used  to  address  the  second  issue,  to  try  to  confirm  the 
mechanical,  essentially  barotropie  nature  of  the  double-gyre 
response.  To  that  end,  the  selected  bora  events  were  studied 
using  barotropie,  wind-only,  single-episode,  QLIODDY 
runs  and  appropriate  segments  of  the  baroclinic,  all- 
fluxes-induded,  long-series,  NCOM  run.  The  output  of 
both  models  was  validated  by  comparison  to  the  JR  P/WISE 
ADCP  current  meter  data.  While  dealing  with  the  third 
issue,  the  offshore  protrusion  of  freshened  l riverine)  and 
cooled  coastal  waters  were  used  to  trace  the  spinup  of  the 
northern  (Trieste)  cyclonic  gyre.  This  part  of  the  study  was 
aided  with  locally  received  AVHRR  data  processed  to 
obtain  analyzed  SST  fields.  The  impact  of  cooling  heat 
exchange  at  the  air-sea  interface  was  also  addressed  in  this 
experiment.  To  that  end,  the  NCOM  16-27  February  2003 
simulated  sequence  of  3-hourly  solutions  provided  the 
reference  case.  The  QUODDY  model,  started  from  homo¬ 
geneous  winter  conditions  (8°C  38  psu)  and  forced  with 
COAMPS  10-m  winds  plus  NCOM -corrected  COAMPS 
surface  heat  fluxes,  provided  the  auxiliary  run  useful  in 
elucidating  the  impact  of  the  heat-flux  inclusion. 

6.  Results  and  Discussion 

6,1.  Variable  Wind  Forcing  and  Gyral  Response 

Pattern 

[2b]  In  the  KO  paper,  the  double-gyre  response  was 
studied  using  a  coarser  grid  and  bathymetry  (7.5-km  reso¬ 
lution),  a  smaller  basin  (down  to  the  Pula-Pesaro  line. 
Figure  1),  and  a  climatological,  steady  wind  forcing.  In 
the  present  study,  a  more  than  twice-larger  modeled  area 
was  covered  with  QUODDY,  and  the  whole  Adriatic  with 
NCOM.  An  order-of-magnitude  finer  grid  was  used  (with 
commensurate  bathymetry)  as  well  as  forcing  with  wind 
akin  to  its  true  variability  in  space  and  time.  The  quality  of 
i  he  COAM  PS  atmospheric  reanalysis  surface  fields,  which 
were  used  to  force  the  sea  models,  is  illustrated  by  com¬ 


parison  of  the  simulated  winds  with  observations  at  three 
meteorological  stations  (Figure  4),  Several  bora-intensity 
oscillations  occurred  during  selected  January  and  February 
episodes. 

[27]  At  the  Pula  station,  the  COAM  PS  simulation  captures 
promptly  and  accurately  the  events  of  the  strengthening 
northeasterly  winds  and  the  relatively  quiescent  periods 
between  events.  The  model  simulations  overall  capture  the 
duration  of  the  strong-wind-speed  event  and  generally  repro¬ 
duce  well  the  magnitude  of  the  events.  However,  the  model 
slightly  under-forecasts  the  magnitude  of  several  of  the  bora 
peaks,  such  as  on  7  and  9  January,  and  over- forecasts  the 
speed  maximum  on  12  January  and  12  February,  The  differ¬ 
ences  between  the  model-simulated  and  observed  maxima 
typically  tall  in  the  range  of  2  5  ms  3, 

[2s]  The  COAMPS  winds  compare  reasonably  well  with 
the  observations  at  Trieste  as  well,  but  the  model  in  general, 
tends  to  overpredict  the  wind  magnitude.  The  model -to-data 
discrepancies  observed  at  several  current  stations,  presented 
later  in  this  section,  support  this  conclusion.  However, 
Signet l  et  ai  [2005]  found  that,  compared  to  w  ind  observed 
at  the  oceanographic  tower  Aqua  Aim  (our  station  VRl), 
COAMPS  predicted  well  a  March  2001  bora,  slightly 
under-estimating  its  magnitude. 

[29]  Notable  di (Terences  between  the  simulated  and 
observed  winds  are  also  observed  at  the  Senj  station. 
For  example,  the  model  predicts  rather  successfully  the 
onset,  duration,  and  magnitude  of  the  9  10  January  bora 
event,  but  then  significantly  over  estimates  the  strength  of 
the  next  peak.  A  similar  effect  is  observed  in  the  February 
episode.  Some  of  these  discrepancies  possibly  arise  from 
small-scale  features  in  the  coastal  topography,  which  are 
likely  important  for  the  local-scale  structure  of  the  bora, 
but  which  may  be  smoothed  or  not  completely  resolved 
in  the  4-k.m  horizontal  resolution  COAMPS.  In  line  with 
that  reasoning  one  notes  that  COAM  PS  excels  at  Pu  la,  which  is 
surrounded  by  low  and  fairly  smooth  orographic  features, 
but  shows  larger  differences  at  Trieste  and  Senj,  where 
orographic  gradients  are  more  pronounced.  Additionally, 
wind  speed  time  series  constructed  for  the  four  COAMPS 
grid  points  surrounding  Senj  (not  shown)  indicate  substantial 
sensitivity. 

[30]  Vorticity  in  the  wind  field  imparts  a  generic  pattern 
on  the  northern  Adriatic  circulation.  Several  requirements 
have  to  be  met  before  this  sea  response  is  properly  simu¬ 
lated.  The  first  and  foremost  is  the  spatial  resolution  of  the 
forcing  w  ind  field.  A  previous  study  [Kuz/nic  et  ai ,  1988] 
provides  evidence  that  it  should  not  be  much  coarser  than 
about  15  km.  Provided  that  the  resolution  requirement  is 
fulfilled,  the  skill  in  predicting  the  onset  and  spatial  extent 
of  the  gyres  and  the  proper  magnitude  and  direction  of  the 
bora-induced  currents  determines  the  eventual  success  of  a 
particular  model  and/or  simulation, 

[31]  It  has  been  shown  in  the  KO  paper  that  the  frictional 
adjustment  times  of  the  shallow,  northern  Adriatic  are  much 
shorter  than  the  typical  duration  of  the  atmospheric  dis¬ 
turbances.  assuring  its  prompt  accommodation  to  the  wind 
forcing.  Furthermore,  lateral  variability  in  the  bora  wind 
field  imparts  a  much  richer  response  than  mere  convergence 
against  the  Italian  coast  would  produce.  A  snapshot  of 
transient  response,  obtained  with  the  COAM  PS-forced 
QUODDY  model,  is  depicted  in  Figure  5  for  both  the 
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Figure  4.  Observed  (solid)  and  CO  AM  PS-predicted  (dashed)  hourly  wind  speed  at  the  Trieste  (top). 
Pula  (middle),  and  Senj  (bottom)  stations  during  the  selected  bora  episodes  in  January  (left)  and  February 
(right)  2003. 


January  (top)  and  February  (bottom)  episodes.  Enhanced 
surface  wind  stress  below  the  Trieste  and  Senj  bora  jets 
generated  southwestward  transports  and  convergence 
against  the  western  coast.  The  convergence  imparts  a  sea- 
level  slope  that  drives  a  southward  How.  However,  the  wind 


vorticity  and  the  related  subsidence  of  the  wind  between  the 
two  jets  (see  the  Rovinj  station  data  in  Figure  2)  modify  the 
convergence.  Consequently,  the  upper  (Trieste),  and  lower 
(Senj)  cyclonic  gyres  are  formed  around  the  corresponding 
lows  of  the  sea-surface  height.  The  basin  east-  and  west-side 
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Figure  5.  A  snapshot  of  the  sea  level  and  circulation  transient  response  to  the  January  (top)  and 
February  (bottom)  bora  episode  as  simulated  with  the  QUODDY  model  started  from  rest  and  forced  with 
COAMPS  winds. 


waters  in  the  area  between  the  two  cyclonic  gyres  respond 
differently  to  the  establishment  of  the  gyres.  In  the  deeper, 
eastern  part,  which  is  bounded  on  its  northern  and  southern 
sides  by  flows  of  opposite  direction,  an  anticyclonic 


(Rovinj)  gyre  is  formed.  In  the  shallower  western  part, 
which  is  under  stronger  slope  control,  and  more  immediate 
Po  River  and  Senj  gyre  influence  more  transient  patterns 
develop.  In  this  area  two  opposing,  southward-  and  west- 
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Figure  6.  Measured  (ADCP)  and  modeled  (NCOM,  QUODDY)  low-pass- filtered  curreni  speed  at 
stations  VRI  (left),  and  VR5  (right)  at  the  depth  of  the  surface  (top),  middle,  and  bottom  cell  location 
during  the  January  bora  episode.  Also  plotted  is  the  COAMPS  wind  stress  at  the  Trieste  station  in  N/nr 
scaled  by  0.5  to  lit  the  current  speed  axis.  The  x  axis  denotes  days  in  January. 


ward-directed  Hows  have  to  accommodate  each  other  (see 
Figure  12  and  discussion  later  on),  and  their  strength  and 
directional  changes  control  the  interplay.  The  Trieste  and 
Senj  cyclonic  gyres  and  the  Rovinj  anticyclonic  gyre  appear 
to  be  p erm a ne nt  tea t u res  o f  t he  respon sc,  with  details  that 


may  vary  from  one  bora  event  to  another.  We  will  examine 
next  from  both  observational  and  modeling  side  the  details 
of  the  response  at  the  ADCP  station  locations,  starting  with 
the  Trieste  gyre. 
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Figure  7.  Progressive  vector  diagram  of  the  vertically  averaged,  low-pass-filtered  currents  at  stations 
VR1  (top  left),  VR2  (top  right),  VR4  (bottom  left)  and  VR5  (bottom  right)  derived  from  the  ADCP 
measurements  (pluses),  the  NCOM  model  (open  circles),  and  the  QUODDY  model  (open  squares)  for  the 
January  bora  episode.  Numbers  mark  the  number  of  days  since  7  January. 


[32]  In  a  shallow,  unstratified  water  column,  wind-stress- 
imparted  momentum  is  transferred  downward  expediently, 
so  currents  along  the  northwestern  coastal  strip  respond 
quickly  to  an  onset  of  bora  wind.  The  simple  analytical 
model  of  Schwing  et  al.  [1985]  for  a  shallow,  coastal  stretch 
of  constant  depth  can  be  used  to  infer  the  phase  relations 
between  the  wind  and  transports.  In  this  frictionally  dom¬ 
inated  situation,  with  an  offshore  distance  of  10  km,  a  depth 
of  1 5  m,  and  a  wind  period  of  2  days,  the  model  predicts 
establishment  of  alongshore  transport  at  VR1  lagging  the 
wind  by  0.5  hrs.  Low-pass-filtered-current  magnitudes, 
observed  and  modeled  at  the  surface  (top),  midcolumn 
(middle),  and  bottom  depths  of  VR1  (left)  and  VR5  (right) 
are  depicted  in  Figure  6.  The  wind  and  currents  appear 
correlated,  and  several  things  are  readily  apparent.  The 
current  response  is  rather  barotropic  and  the  bora  signal  is 


transmitted  down  to  the  bottom.  The  models  reproduce  that 
with  remarkable  success  regarding  the  timing  and  magni¬ 
tude  of  the  bora-induced  peaks,  relative  to  the  observations. 
The  two  models  forced  with  the  same  wind,  but  otherwise 
different  in  many  design  and  forcing  aspects,  produce 
current  magnitude  predictions  generally  following  each 
other.  Both  models  underpredict  the  9  10  January  peak, 
and  drastically  overpredict  the  1 1  13  January  peak.  One 
notes  at  VR5  better  correspondence  between  the  observed 
and  modeled  peak  magnitudes  and  timing,  for  the  period 
11-13  January  in  particular.  Progressive-vector-diagram 
(PVD)  plots  of  the  vertically  averaged  currents  at  VRI 
and  VR2  for  the  same  episode  (Figure  7,  top)  offer 
additional  insight.  The  plots  provide  pseudo-Lagrangian 
information,  which  is  more  susceptible  to  directional  differ¬ 
ences.  In  accord  with  the  shallowness  and  lateral  con- 
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straints,  the  measured  and  modeled  paths  at  VRI  and  VR2 
closely  follow  each  other  throughout  the  six-day  period  (the 
ADCP  and  QUGDDY  path  at  VR2  in  particular).  The 
longer  modeled  paths  at  both  locations  are  a  consequence 
of  over-predicting  the  current  speeds.  The  February  episode 
(not  shown)  pictures  essentially  the  same  situation  with  20 
30  km  shorter  paths. 

[33]  The  old  KO  simulations  showed  that  in  a  flat-bottom 
non  hem  Adriatic,  Coriolis-deflected  sea-surface  slope  is 
established  within  a  day*  When  the  basin  is  driven  by  a 
horizontally  varying  wind  field,  variations  in  both  depth  and 
wind  stress  become  important.  Bathymetric  and  wind 
effects  can  be  evaluated  using  the  linearized,  conservation 
of  mean  flow  vorticity  equation  [Csatwdv,  1975]: 


where  u  and  v  are  depth -averaged  velocities,  r  is  the  wind 
stress  vector,  and  h  is  the  water  depth.  If  the  relati  ve  change 
of  the  wind  stress  is  greater  than  the  topographic  variation 
(which  is  the  case  in  the  northern  Adriatic),  the  wind 
variations  dominate.  As  pointed  out  a  long  time  ago 
[Shtokman,  1941],  in  such  situations,  lateral  heterogeneities 
in  the  wind  field  produce  circulations  that  close  both 
horizontally  and  vertically,  resulting  in  an  upwind  surface 
flow  at  some  locations.  Such  a  situation  is  observed  at 
stations  VR4  and  VR5.  In  contrast  to  VRI  and  VR2,  where 
the  slope  currents,  developed  in  response  to  the  cross-shore 
Coriolis  force,  reinforce  the  drift  currents,  at  VR4  and  VR5 
the  slope  currents  oppose  a  weaker,  stress- imparted  drift  and 
the  net  flow  is  upwind  (Figure  7,  bottom).  One  readily  notes 
that  at  both  VR4  and  VR5  model  predictions  closely  follow 
the  observations  until  10  January,  and  disagree  more  in  the 
second  part  of  the  January  episode.  Inspection  of  Figure  4 
shows  generally  better  COAMPS-magnitude-to-data  agree¬ 
ment  during  the  first  part  of  the  January  episode, 

[34]  The  area  southeast  of  the  Trieste  gyre  responds  to  the 
bora  with  a  cell  circulating  in  the  opposite  direction  (see 
Figure  5),  But  here  the  intensity  of  the  currents  is  much 
lower  and  the  model/observations  discrepancy  is  much 
higher  than  at  the  stations  under  direct  bora  forcing.  At 
both  the  VR6  and  IC  stations  and  for  both  bora  episodes 
(with  the  partial  exception  of  VR6  in  January),  the  current 
data  dearly  trace  a  part  of  the  southeasterly  directed 
(antieyckmic)  flow  (Figure  8).  In  terms  of  the  vertically 
averaged,  low-pass- filtered,  progressive  current  vectors,  the 
model-to-data  agreement  is  highest  at  VR6  in  February  and 
lowest  at  IC  in  January.  The  agreement  between  the  models 
(in  spile  of  the  same  wind- field  forcing)  is  much  lower  than 
at  the  other  VR  stations,  with  occasional,  completely 
opposite  direction  of  the  predicted  flow  (at  IC  in  January). 
Very  similar  PVD  paths  are  obtained  for  the  surface  bin 
positions  at  both  stations  and  for  both  episodes  (not  shown), 

[35]  Modeled  and  observed  speed  at  CP3  and  CP5  during 
the  January  episode  is  plotted  in  Figure  9.  As  in  Figure  6, 
the  surface  response  is  plotted  at  the  top.  the  midcolumn 
response  in  the  middle,  and  the  bottom-cell  level  response 
speed  in  the  lowest  row.  The  CP3  station  exhibits  currents 
of  lower  intensity  than  at  VRI,  but  both  models  as  before 
predict  mutually  similar  speed  values.  The  most  problematic 


again  is  the  end  of  January  episode,  with  both  models 
notably  overpredicting  the  speed.  The  situation  is  worse  at 
CPS  where  both  models  first  underpredict  the  speed  and 
then  grossly  overpredict  it  by  the  end  of  the  episode. 
Inspection  of  the  Senj  wind  in  Figure  4  points  to  the  wind 
field  as  the  cause  of  the  discrepancy.  The  PVD  paths  in 
Figure  10  for  CP3  (middle),  and  CPS  (bottom)  show  the 
models  overpredicting  the  northward  How  at  CPS  and  the 
southward  flowatCP3.  At  the  westernmost  station,  CP2,  both 
models  predict  a  steady  southward  movement  at  odds  with  the 
ADCP  data,  which  exhibits  much  slower  eastward  flow. 

[3ft]  In  order  to  gain  additional  insight  and  to  compare 
more  formally  the  differences  between  the  measured  and 
modeled  response,  several  useful  ratios  are  introduced 
[Schwab  er  aL  I989J.  Specifically: 

-Ik  -  u*):+(v*  -  v*r] 


where  rx  is  the  ratio  of  the  time-variant  current  energy  to  the 
total  current  energy  with  x  equal  to  O(bserved),  N(COM)  or 
Q(UODDY)  summation  is  over  hourly  values  within 
selected  period: 


E 

[[(«*  -«»)  -•  («v  -  Mf+tK  ■ 

-  (v>  -  v>)]'j 

f*V\y  — 

=1 

[cu«  -  ■u«)I+(v,  -  v,r 

1 

(6) 


(U\  -  Uv}  +(VX  -  Vv 


(?) 


where  rxv  is  the  ratio  of  the  root  mean  square  (RMS)  of  the 
computed  currents  to  the  RMS  of  the  observ  ed  currents. 
is  the  ratio  of  the  computed  and  observed  total  flow.  n-u  is 
the  ratio  of  the  computed  and  observed  time- variable  flow, 
and  nnxy  is  the  ratio  of  the  computed  and  observed  time- 
mean  flow;  the  xy  equals  ON  or  OQ.  The  comparison  is 
summarized  in  Tables  3a  -  3c.  At  VRI,  both  the  observed 
and  modeled  rK  values  are  less  than  0,2  for  both  episodes 
(with  the  exception  of  r<>  for  January),  indicating  high 
values  of  episode-averaged  mean  energy.  At  CP5,  the  upper 
limit  for  ?\  is  almost  twice  as  high  (0.388),  suggesting  a 
proportionally  lower  episode-averaged  mean.  The  t\  values 
are  highest  at  IC  where  episode  and  depth-averaged  speeds 
exhibit  the  lowest  values  and  time-variable  energy  can 
account  for  more  than  80%  of  the  total  energy.  Both  models 
underpredict  this  ratio  for  both  episodes  at  the  directly 
wind-forced  stations  (VRI  and  CP5),  with  the  exception  of 
QUODDY  at  CP5.  At  1C,  the  models  exhibit  n  and  rv  ratios 
greater  than  one,  with  NCOM  providing  a  factor  of  3  5 
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Figure  8.  Progressive  vector  diagram  of  the  vertically  averaged,  low-pass- filtered  currents  at  stations 
VR6  (left)  and  IC  (right)  for  the  January  (top)  and  February  (bottom)  episodes.  Numbers  mark  the 
number  of  days  since  7  January  and  9  February,  respectively. 


larger  variability  than  observed  in  the  data.  The  RMS  values 
of  the  computed  currents  are  always  higher  than  the 
respective  observed  values  for  NCOM  (indicating  over 
prediction)  and  lower  than  observed  tor  QUODDY  (with 
the  exception  of  VR I  for  January  and  1C  for  February). 

[u J  Calculated  observed/model  differences  for  the  total, 
time- variable,  and  mean  flow  provide  additional  insight.  At 
all  stations,  QUODDY  gives  smaller  rt  values  than  NCOM 
for  both  episodes,  suggesting  smaller  depth-  and  episode- 
averaged  discrepancies  with  the  observed  total  flow.  The 
same  is  true  for  both  the  time- variable  and  mean  component 
of  the  flow,  writh  the  exception  of  the  VR  1  value  for  January. 
Even  that  exception  appears  to  be  related  to  anomalously 
low  rv0N  values  near  bottom.  The  QUODDY-calculated 
mean  currents  agree  rather  well  with  the  observed  means  at 
VRI  and  CP5  and  clearly  not  so  well  at  1C.  At  all  three 
stations  and  for  both  episodes,  the  QUODDY  differences 


with  the  observed  are  variously  smaller  than  NCOM.  One 
possible  explanation  is  that  NCOM  accumulated  over  the 
previous  simulated  months  a  component  of  the  flow,  for 
wfhich  the  episode-focused  QUODDY  did  not  have  the  time. 
Discrepancies  in  the  calculated,  time-variable  pans  of  the 
currents  exceed  the  respective  observed  variability  for  both 
models  and  all  stations  and  episodes;  for  NCOM  the  rvx> 
ratio  reaches  a  value  of  5.70  at  the  1C  surface  in  the  January 
episode  and  drops  down  to  1.20  at  the  VRI  bottom  in  the 
February  episode.  The  QUODDY  values  are  always  lower 
and,  apan  from  the  stated  extremes,  follow  NCOM.  When 
discussing  these  values,  one  has  to  bear  in  mind  that 
observed  values  as  processed  may  still  contain  current 
fluctuations  beyond  the  models'  predictive  skills.  We  did 
calculate  the  same  statistics  with  an  Sih-order,  Butterworth, 
low-pass  filter  with  a  cutoff  period  of  36  hours.  Filtering 
with  a  longer  period  predictably  lowered  the  r(>  values  by 
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Figure  9,  Measured  (ADCP)  and  modeled  (NCOM,  QUODDY)  low-pass-filtered  current  speed  at 
stations  CP3  (left)  and  CP5  (right)  at  the  depth  of  the  surface  (top),  middle,  and  bottom  cell  location 
during  the  January  bora  episode.  Also  plotted  is  the  COAMPS  wind  stress  at  the  Trieste  station  in  N/rrr 
scaled  by  0.5  to  fit  the  y  axis.  The  x  axis  denotes  days  in  January. 


a  few  percent,  but  other  statistics  did  not  improve*  Also,  it 
is  worth  reiterating  here  that  the  NCOM  simulation  is  a 
long*  seven-month  run  with  momentum*  surface  heat*  and 
tidal  forcing  (filtered  out  before  comparison),  whereas 


QUODDY  provides  only  episodic  simulations  of  the  wind- 
forced,  barotropic  sea.  The  obtained  metric  differences 
invoke  discussion  of  deeper  issues  beyond  the  scope  of  this 
paper,  like  deterministic/stochastic  nature  of  the  modeled/ 
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Figure  10,  Progressive  vector  diagrams  of  the  vertically 
averaged,  low-pass-filtered  currents  at  stations  CP2,  CPS, 
and  CP5  derived  from  the  ADCP  measurements  (pluses), 
NCOM  model  (open  circles),  and  QUODDY  model  (open 
squares)  for  the  January  bora  episode.  Numbers  denote  the 
number  of  days  since  7  January, 

observed  variability,  or  high-resolution  model  errors  at  low 
frequency  and/or  low  wave  number  scales.  In  their  study  of 
wind  quality  for  oceanographic  modeling.  Signet!  et  at 
[2005]  derived  surface  winds  from  four  meteorological 
models  obtaining  better  model- to-data  correlation  with 
coarser  model  winds.  In  another  study  quoted  in  their  paper 
[Bogden  et  at 1996]  a  2D  wind-forced  linear  barotropic 


model  outperformed  a  more  realistically  forced  3D 
nonlinear  model  based  on  a  mod  el /data  misfit  metric. 
One  is  well  advised  not  to  take  the  obtained  NCOM- 
to-QUODDY  differences  to  represent  universal  skills  of 
the  models  in  question. 

[Jfi]  Akin  to  the  data-only  analysis  presented  in  Figure  3, 
the  north-south  components  of  both  the  measured  and 
modeled,  vertically  integrated  and  low -pass- tillered  veloci¬ 
ties  at  VRL  1C,  and  CP5  are  plotted  as  scatter  diagrams 
against  their  east-west  counterparts  in  Figure  11.  At  VRL 
the  models  for  both  episodes  predict  very  well  a  downwind, 
coast-constrained  flow.  The  polarization  is  sharper  than 
observed  and  average  speeds  are  over  predicted.  At  CP5, 
the  data  exhibit  less  polarization  and  the  models  to  some 
degree  capture  that  change.  The  weak  response  at  station  JC 
has  proved  more  difficult  to  reproduce.  In  particular,  both 
models  mispredict  the  orientation  of  the  principal  variability 
axis  by  several  tens  of  degrees.  The  level  of  the  models'1 
skill  in  predicting  the  northern  Adriatic  response  to  strong 
bora  events  can  be  still  better  appreciated  by  looking  at  the 
depth-  and  episode-averaged  flows  (Figure  12).  The  mod¬ 
els'  responses  as  well  as  measurements  plotted  in  the  figure 
were  first  averaged  over  the  time  span  of  the  selected 
episodes,  and  then  over  the  local  depth.  The  frictional 
adjustment  time  being  much  shorter  than  the  atmospheric 
disturbance  (KG)  assures  a  swift  response  and  justifies  the 
averaging  period.  Clearly,  both  models  and  both  bora  events 
produce  the  same  basic  gyral  pattern.  The  observed  and 
modeled  flows  exhibit  high  similarity  in  position,  magni¬ 
tude  and  direction,  but  there  are  some  differences.  Judging 
by  the  VR5  data  (see  Figure  I )  both  models  appear  to 
underpredict  the  upwind  flow,  in  February  episode  in 
particular.  For  the  January  episode  both  models  reproduce 
somewhat  less  convincingly  the  flow  direction  at  the  two 
stations  covering  the  antteydonic  gyre  (VR6  and  1C)  than  at 
VR1  and  VR2,  but  nevertheless  model-derived  averaged 
velocities  compare  well  with  equally  averaged  measured 
values,  fn  all  the  four  averaged  fields  one  notes  an  area  of 
sluggish  currents  in  which  two  opposing  and  bifurcating 
flows  meet:  one  heading  south,  reinforcing  the  Western 
Adriatic  Current  (WAC),  and  the  other  moving  westward, 
branching  off  from  the  Senj  gyre  northern  arm.  The  location 
of  this  area  is  somewhat  different  in  both  models,  and  for 
both  episodes,  always  positioned  nonh/nonhwest  of  the 
CP2  station,  where  measurements  indicate  very  small  aver¬ 
age  current. 

[39]  To  further  evaluate  the  model/data  discrepancies  in 
the  gyral  response,  the  magnitude  of  the  complex  correla¬ 
tion  coefficient  [p)  and  the  angular  displacement,  or  mean 
directional  error,  (<p)  between  the  measured  (ADCP)  and 
modeled  (NCOM.  QUODDY)  currents  were  calculated 
following  Kundu  [1976]: 


(tipUm  +  \>vm)  T  f{u$vm  umv0) 

\J (ur>  +  Vn  )  ■  \J  { um  +  vm  ) 


tp  —  arerg 


{u(>vni  umvti) 

<u„um  +  V0Vffl) 
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where  n  and  v  are.  respectively,  east-west  and  north-south, 
observed  {o>  or  modeled  (m),  demeaned  velocity  compo- 
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Table  3a.  Comparative  Low-Pass-Filtered  Current  Statistics:  VRU 


Period 

Depth 

To 

tn 

ro 

0>N 

Too 

rtoN 

rtoo 

rv()N 

r vtK> 

nn„N 

nniKj 

Jan 

1.7 

.236 

189 

190 

1.856 

1 .540 

1.222 

,778 

2.473 

2.182 

,835 

.346 

episode 

6.2 

.243 

.177 

. !  89 

1.621 

1.439 

0.849 

,693 

1,994 

2.096 

.482 

244 

14.7 

.252 

.157 

.191 

1.350 

1.369 

0.487 

637 

1.350 

1,951 

.195 

.193 

Vbar 

.244 

174 

,190 

1,574 

1,4 1 7 

.799 

.677 

1.929 

2,072 

.438 

,228 

Feb 

1.7 

.086 

.098 

.141 

1  255 

.915 

.308 

.215 

2,926 

2.337 

.063 

.016 

episode 

6.2 

.1 14 

.091 

.138 

1  145 

.879 

.273 

.235 

2.167 

1,889 

.028 

.022 

14.7 

.227 

.082 

.135 

1.002 

.875 

.286 

.275 

1,254 

1.188 

.1108 

.007 

Vbar 

.145 

.090 

.137 

1.1 15 

.872 

.279 

.247 

1.975 

1.691 

.028 

,02 1 

J  FI  ere  r( 

is  the  ratio  of  the  time-variant  current  energy  to  the  lotal  current  energy  with  x  equal  to  Observed),  N(COM)  or  QiUODDY);  r„  is  the  ratio  of 

ihc  root  mean  square  (RMS)  of  the  computed  currents  to  the  RMS  of  the  observed  currents.  rty)  is  the  ratio  of  the  computed  and  observed  total  flow,  n\L  is 
the  ratio  of  the  computed  and  Observed  time-variable  flow,  rm „  is  the  ratio  of  the  computed  and  observed  time-mean  flow:  the  xy  equals  ON  or  00;  vbar 
marks  the  column  averaged  values.  The  ratio  equations  are  provided  in  the  text;  summation  is  over  hourly  values  within  selected  period.  Depth  is  in  meters 


nents  and  the  brackets  represent  time  averaging.  For  each 
model  and  both  time  periods,  the  correlations  were 
computed  tor  each  station  as  a  function  of  depth,  starting 
at  the  topmost  ADCP  bin.  The  complex  correlations 
between  the  ADCP -observed  and  NCGM/QUODDY  model- 
predicted  currents  were  typically  rather  uniform  throughout 
the  water  column.  As  a  next  step,  both  p  and  (p  were  further 
depth  averaged  to  provide  one  pair  of  indicators  per  station 
and  episode.  The  amplitude  correlation  results  are  shown  in 
Figure  13.  The  best  amplitude  correlation  was  obtained  tor 
VR5  in  January  and  the  worst  at  VR1  in  February,  It  is 
clearly  related  to  the  fact  that,  during  the  January  episode, 
both  model  predictions  closely  follow  the  measured  speed  at 
VR5  (Figure  6),  whereas,  during  the  February  episode  at 
VR1  (not  shown),  there  was  a  continuous  mismatch  in 
speed,  with  the  models  first  over  predicting  the  observed 
values  and  then  under  predicting  them  by  the  end  of  the 
episode.  As  expected  from  the  presented  PVDs  and 
scatterplots,  both  models  generated  small  (10°)  average 
veering  at  VR 1  for  the  January  episode,  but  larger  directional 
discrepancies  at  CP5- 

6,2.  Bora-Driven  Extrusion  of  Coastal  Waters 

[40]  The  third  issue  we  wish  to  explore  can  be  traced  back 
to  the  Kuzmlc  [1991]  paper  in  which  virtual  particles  and 
reflectance-derived  chlorophyll  concentration  were  used  to 
map  the  horizontal  extent  of  bora- induced  return  flow'  in  a 
simple  kinematic  framework.  Alternatively,  the  SST  can  be 
used  to  trace  bora-induced  flow  features.  We  have  confined 
simulations  reported  in  this  section  to  the  case  of  cooling  as 
a  source  of  cold  water  usable  as  that  kind  of  tracer. 

[41 ]  A  flux  coupler  option  was  implemented  in  NCOM  to 
calculate  latent  and  sensible  heat  fluxes  using  bulk  formulas 


with  NCOM-predicted  SST  and  COAM PS-derived  atmo¬ 
spheric  fields.  The  very  same  heat  flux  was  also  used  in  a 
related  simulation  to  drive  the  QUODDY  model 

[42]  This  experiment  is  focused  on  the  16-27  February 
time  period  when  AVHRR  scenes  were  available.  Develop¬ 
ment  during  this  twelve-day  sequence  is  preconditioned  in 
the  sense  that  it  follows  the  strong  February  bora  episode 
studied  in  the  previous  section,  in  those  analyses,  however, 
the  thermal  aspect  was  not  considered.  To  address  this,  we 
ran  an  additional  QUODDY  simulation  covering  the  period 
8  16  February  since  for  these  terminal  dates,  relatively 
clear  AVHRR  scenes  were  available  for  comparison.  More 
specifically,  to  aid  the  NCOM  output  analysis,  the 
QUODDY  model  was  forced,  tor  the  same  period,  with 
the  same  NCOM -corrected  heat  flux  starting  from  a  homo¬ 
geneous,  ditnatolpgically  correct,  initial  state  (T  ■=  8°C.  S  = 
38  psu).  The  outcome  is  summarized  in  Figure  14,  The  SSI 
difference  (16  February-8  February)  as  obtained  wuih 
the  NCOM  model  is  given  in  the  bottom  left  corner 
(Figure  14c),  and  the  same  output  from  the  QUODDY 
model  is  plotted  in  the  bottom  right  comer  (Figure  14d), 
The  two  mentioned  AVHRR  scenes  are  also  included  for 
reference;  one  of  which  covers  the  beginning  (Figure  14a), 
and  the  other  the  end  (Figure  14b)  of  the  simulated  period. 
Due  to  considerable  cloudiness  on  both  images  (which  adds 
up  in  the  result)  differencing  of  the  images  was  done,  but  is 
not  presented.  Unobstructed  parts  of  the  AVHRR  difference 
image  suggest  a  similar  situation  to  that  observed  in  the 
NCOM  and  QUODDY  difference  plots,  he.,  a  lowering  of 
the  average  northern  Adriatic  temperature,  a  colder  strip 
along  the  shallow,  northwestern  coast,  and  an  extruded 
tongue  of  colder  water.  Whereas  the  NCOM  difference  is 
closer  to  reality  (e.g„  exhibiting  islands  of  wanner  water 


Table  3b,  Comparative  Low- Pass- Filtered  Current  Statistics:  1Ca 


Period 

Depth 

To 

r,M 

Tfi 

TOO 

ttn\ 

rtoo 

rVQN 

rvt*> 

mio* 

Jan 

7.2 

.264 

.669 

Ml 

1 ,683 

.788 

4.706 

1.108 

5.647 

0.902 

4,351 

1.1X2 

episode 

19.2 

.304 

.790 

.838 

1.280 

,688 

3.063 

0.830 

3.440 

0.785 

2.898 

0.850 

39.2 

.284 

.947 

.820 

1.155 

.750 

1,765 

0.700 

4.43N 

1.379 

0,497 

0.428 

vbar 

.294 

.853 

.816 

1.259 

,717 

2.849 

0.818 

4.145 

0.986 

2.320 

0.751 

Feb 

7.2 

.690 

4X1 

,424 

1 ,614 

1,473 

3.484 

2.305 

3.271 

2.097 

3,956 

2.767 

episode 

19.2 

.762 

.600 

*565 

1458 

1.240 

2,493 

2.134 

2.254 

l  ,944 

3,257 

2.744 

39.2 

.560 

*595 

.673 

1.010 

1.044 

1.246 

1,665 

1 ,5K6 

2.313 

0.813 

0.839 

vbar 

.694 

*566 

,589 

1.167 

1.202 

2.265 

2.007 

2.136 

2.059 

2-703 

2.109 

3 As  in  Tabic  3a. 
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Comparative  Low-Pass-Filtered  Current  Statistics:  CP5“ 

C03SI3 

Period 

Depth 

ro 

Tn 

ro 

Ton 

TfJO 

^ON 

rtoo 

r Vqn 

rv \  )U 

nnON 

nil* 

Jan 

7.9 

.294 

.259 

305 

1.399 

.983 

1.271 

.*■13 

3,651 

2.666 

378 

.041 

episode 

21.9 

.316 

.273 

344 

1.233 

.898 

1.099 

,800 

3. Ml 

2.421 

170 

.052 

39.9 

339 

.334 

381 

1.003 

.804 

3.038 

.768 

2,832 

2.080 

119 

.096 

vbar 

.322 

385 

.349 

1.199 

.879 

1.100 

.787 

3.100 

2336 

167 

.062 

Feb 

7,9 

.223 

339 

.316 

1.854 

1  155 

1294 

.372 

2.735 

1.476 

,752 

.055 

episode 

21.9 

332 

.276 

350 

t.628 

L010 

0.801 

.329 

2.200 

1,264 

379 

U47 

39.9 

33  8 

.388 

.380 

1370 

0.848 

0.491 

392 

1.853 

1.269 

.066 

.119 

vbar 

.236 

.295 

.354 

1.576 

0.983 

0.774 

350 

2,144 

1369 

354 

.067 

JAs  in  Table  3a, 


advected  from  the  south)  the  QUODDY  difference  is  Po  River  was  not  included  in  the  QUODDY  simulation 
perhaps  more  instructive.  The  2.5  -  J.0°C  temperature  drop  either  as  a  flow/velocity  or  temperature  boundary  condition, 
observed  in  the  plot  clearly  traces  the  extent  of  both  the  Therefore,  the  tongue  of  extruded  colder  water  seen  in  the 
Trieste  and  Senj  gyres.  It  is  important  to  stress  here  that  the  QUODDY  differential  image  is  the  sea  water  cooled  along 
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Figure  II.  Scaiterplot  of  the  vertically  averaged  currents  at  the  VRI,  1C,  and  CPS  stations:  modeled 
with  NCOM  (left),  ADCP  measured  (middle ),  and  modeled  with  QUODDY  (right)  during  the  January 
(top)  and  February  (bottom)  bora  episodes. 
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NCOM  QUODDY 


Figure  12.  Depth  and  episode  averaged  flow  computed  with  the  two  models  and  measured  at  nine 
stations.  Thick  vectors,  with  dots  at  their  roots,  mark  the  averaged  measured  currents.  The  NCOM  output 
is  rotated  to  the  northeast  grid  and  resampled  with  0,1°  resolution.  The  unstructured -mesh  QUODDY 
output  is  mapped  onto  the  same  grid.  The  time  averaging  periods  equal  the  duration  of  the  selected  bora 
episodes. 


the  shallow,  northwestern,  coastal  strip  and  spun  up  offshore 
by  the  bora  wind.  Another  cold-water  tongue  is  obtained 
further  south  tracing  the  northern  arm  of  the  Senj  gyre. 
This  cooling,  traceable  as  well  in  the  AVHRR  scenes  in 
Figure  14,  entails  intensification  of  the  meridional  gradient 
in  the  area  south  of  the  Istrian  peninsula,  a  result  also 
obtained  by  B,  Cush  man -Roi  sin  and  K.  A.  Korotenko 
(Mesoscale-resolving  simulations  of  summer  and  winter 
bora  events  in  the  Adriatic  Sea,  submitted  to  Journal  of 
Geophysical  Research,  2006)  in  their  mesoscale-resolving 
simulation  of  the  winter  bora  (11  20  February  2003). 

[4.']  The  16-27  February  simulation  was  done  only  with 
the  NCOM  model  (because  it  included  the  thermal  fluxes  at 
the  air-sea  interface).  The  AVHRR  data  and  this  simulation 
outcome  are  presented  in  Figures  15  and  16,  respectively. 


The  NCOM  simulation  of  the  northern  Adriatic  SST  over 
the  twelve-day  period  (Figure  16)  is  to  be  compared  to  the 
same  sequence  derived  from  the  AVHRR  observations 
(Figure  15).  The  same  palette  is  applied  in  both  sequences. 
On  top  of  each  figure,  the  wind  magnitude  measured  at 
Trieste  station  is  plotted  to  aid  the  analysts.  The  stars  on  the 
plot  mark  the  slightly  variable,  midday  times  when  the 
AVHRR  scenes  were  acquired  and  the  vertical  bars  mark 
the  noon  of  each  day  for  which  NCOM  predictions  were 
calculated.  Although  the  AVHRR  scenes  are  contaminated 
with  clouds,  both  sequences  clearly  exhibit  a  meridional 
temperature  difference  and  a  frontal  zone  approximately 
following  the  50-m  bathymetry  line.  This  frontal  zone 
separates  the  shallower  north,  which  is  impacted  by  fresh¬ 
water  runoff  and  more  pronounced  cooling,  from  the  deeper 
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Figure  13,  Vertically  averaged  magnitudes  of  the  complex  correlation  between  the  low-pass-filtered 
model-predicted  (NCOM  QUODDY)  and  observed  currents  at  nine  stations  during  the  January  (top)  and 
February  (bottom)  bora  episodes. 


south,  which  is  influenced  by  adveetion  of  warmer  and 
saltier  waters.  Zore- Arman  da  and  Gacic  [19R7]  related  the 
position  of  this  thermal  front  to  the  frequency  of  the  bora 
wind. 


[44]  In  Figure  15  one  notes  on  16  and  17  February 
accumulation  of  colder  water  south  of  the  Po  delta;  the 
water  is  both  discharged  by  the  Po  River  and  adveeted  there 
alongshore  by  the  Trieste  gyre.  The  bulge  appears  colder 
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Figure  14*  Sea  surface  temperature  observed  with  the  AVHRR  sensor  on  (a)  8  February  2003  and 
(b)  16  February  2003.  The  plots  al  the  bottom  show  the  related  sea  surface  temperature  difference 
(16*02,2003-08.02.2003)  as  calculated  with  the  (c)  NCOM  and  (d)  QUODDY  models. 


(and  closer  to  the  4.2°C  and  3.5°C  temperature  values 
measured  on  16  and  17  February,  respectively,  at  a  Po 
upstream  location  -J.  Chiggiato.  personal  communication) 
in  the  AVHRR  scenes  than  in  the  NCOM  simulations*  The 
difference  should  be  partly  ascribed  to  the  fact  that  no  actual 


river  temperatures  were  used  as  a  Po  River  boundary 
condition  in  the  NCOM  simulation  (Adriatic  monthly 
climatological  values  were  used  instead).  Consequently, 
one  can  see  on  most  of  the  images  in  Figure  16  tiny  islands 
of  anomalously  higher-temperature  water  at  ihe  locations  of 
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Figure  15,  Sea  surface  temperature  during  the  period  16-27  February  2003  as  observed  with  the 
AVHRR  sensor  Plotted  at  the  top  is  the  wind  magnitude  measured  at  the  Trieste  station.  Stars  mark  the 
times  of  the  AVHRR  scenes,  and  vertical  bars  mark  the  times  of  the  corresponding  SST  fields  as 
predicted  with  the  NCOM  model 
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Figure  16.  Sea  surface  temperature  during  the  period  16-27  February  2003  as  predicted  with  the 
NCOM  mode).  Plotted  at  the  top  is  the  wind  magnitude  measured  at  the  Trieste  station.  Vertical  bars  mark 
the  times  of  the  model-predicted  SST  fields,  and  stars  mark  the  registration  times  of  the  corresponding 
AVHRR  scenes, 
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the  Po  River  outtlow.  One  should  bear  in  mind  that  this 
sequence  also  exhibits  the  remnants  of  the  previously  spun 
up  colder  water  visible  in  both  the  remotely  sensed  image 
and  the  NCOM  simulation  for  16  February.  The  model 
appears  to  predict  a  sharper  than  observed  frontal  gradient 
around  9°C.  Although  clouds  obstruct  the  view  of  the 
northwestern  shallow  strip  throughout  most  of  the  sequence, 
closer  inspection  reveals  fairly  well  the  correspondence 
between  the  modeled  and  observed  SST.  On  the  background 
of  the  previous  thennal  signature,  one  observes  the  extru¬ 
sion  of  new  cold  water  spun  up  by  the  two-day  bora  event, 
w  hich,  after  a  break  of  about  a  day  and  a  half,  followed  our 
February  episode.  The  extrusion  can  be  traced  in  both  the 
observ  ed  and  simulated  sequences.  Kuzmic  [1991]  exploited 
a  circumstance  that  these  waters  under  riverine  influence 
(the  Po  River  in  particular)  are  also  highly  turbid  and 
therefore  also  observable  by  optical  sensors  (CZCS). 
Present-day  sensors  are  even  more  useable  in  that  regard, 
as  the  Figure  1.4  in  Lee  et  al.  [2005b],  for  example,  clearly 
testifies.  The  reported  SeaWiFS  image  registered  on 
23  February  2003  maps  out  a  chlorophyll  concentration 
pattern  very  similar  to  the  SST  pattern  presented  in  Figure  1 5. 
With  the  bora  subsiding  over  the  following  days,  subsequent 
images  of  the  cold  water  tongue  exhibit  instabilities  and 
meandering  (more  pronounced  in  the  AVHRR  images  than 
in  the  NCOM  simulation)  and  eventual  dissolution  of  the 
temperature  tongue.  Throughout  this  period,  one  can  trace  in 
the  frontal  zone  on  the  southern  boundary  of  the  basin  the 
creation  and  destruction  of  wann  water  filaments.  A  series 
observed  over  a  week-long  period  (21  -27  February)  in  the 
AVHRR  scenes  appears  to  trace  development  of  a  filament 
controlled  by  a  clockwise  motion.  The  same  series  is  visible  but 
less  pronounced  in  the  NCOM  sequence,  which,  over  a 
somewhat  shorter  period  (24-27  February),  shows  the 
development  of  a  similar  filament  closer  to  the  tip  of  Istria 
and  the  Rovinj  anticyclonic  gyre. 

7.  Summary  and  Conclusions 

[4^]  A  combination  of  recent  intensive  observations  and 
simulations  with  two  numerical  models  was  used  to  revisit 
the  issue  of  the  northern  Adriatic  response  to  strong  bora 
episodes.  The  simplistic  simulations  reported  in  the  KO 
paper  have  been  reconsidered  in  the  light  of  recent  field 
evidence  and  more  insightful  modeling  results.  This  new 
information  is  summarized  in  Figure  17b.  in  a  manner 
convenient  for  synthesis,  and  comparison  with  previous 
results.  Both  new  observed  and  simulated  data  reinforce 
the  KO  finding  that  an  episode  of  strong  bora  wind 
provokes  a  double-gyre  (cyclonic.  Trieste,  and  anticyclonic. 
Rovinj)  response  of  the  northern  Adriatic  north  of  the  Po 
Delta  -  Pula  line  (Figure  17b).  The  old  KO  simulations, 
based  on  a  linear  model  and  steady  climatological  wind 
forcing,  successfully  predicted  the  mechanistic,  barotropic 
nature  of  the  induced  gyral  pattern  and  the  appearance  of  an 
upwind  return  flow  between  the  counter-rotating  gyres 
(Figure  17a).  New  modeling  solutions,  driven  by  higher- 
resolution.  time-variant  wind,  position  the  Trieste  cyclonic 
gyre  predominantly  north  of  the  Po  Delta  Rov  inj  line,  in 
good  agreement  with  new'  empirical  evidence  along  the 
gyre's  northeastern  (NE)  and  northestem  (NW)  arms  (area  1 
in  Figure  17b).  The  gyre's  NW  ami  is  geometrically  con- 
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Figure  17.  Schematic  of  the  bora-induced  circulation  in 
the  northern  Adriatic:  (a)  as  reported  in  the  KO  paper;  (b)  as 
derived  in  the  present  work.  The  numbers  in  Figure  17b 
mark  the  four  characteristic  bora-response  areas  discussed 
in  the  text.  Also  plotted  in  Figure  17b  are  the  depth-  and 
episode-averaged  currents;  thick  vectors,  January  episode; 
thin  vectors,  February  episode. 

strained  by  the  northwestern  Adriatic  coast  and  intensified 
by  the  shallowness  of  the  coastal  zone.  During  an  intense 
bora  episode,  both  the  measured  and  modeled  statistics 
picture  a  downwind,  highly  polarized  (£  <  0. IS),  and 
almost  depth-independent  (*>  >  0.93)  flow  within  the  ami. 
The  other  (NE)  arm  maintains  a  sharp  polarization  and  high 
depth  dependence  while  exhibiting  lower  speed,  with  both 
models  again  in  good  accord  with  observations.  Although 
there  were  no  ADCP  measurements  in  the  gyre's  southern 
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section,  water  movements  deducible  from  a  sequence  of 
AVHRR  SST  scenes  reinforce  the  model  solutions. 

[46]  Both  the  simulations  and  measurements  place  the 
other.  Rovinj  gyre  (area  2  in  Figure  17b)  off  the  west  Istrian 
coast,  with  no  comparable  field  information  across  the 
basin,  and  closer  to  the  other  shore.  The  current  statistics 
for  the  Rovinj  gyre  provide  lower  maximum  and  average- 
speed  values  and  less  polarized  (f  <  0.54)  but  still  rather 
depth  independent  (y  >  0.94)  flow,  while  predicting  a 
clockwise  rotation.  Apart  from  some  directional  discrep¬ 
ancies,  the  gyre  is  clearly  in  accord  with  measurements. 
South  of  the  Po  Delta  -  Pula  line  (area  3),  both  models 
predict  a  southward  flow  contributing  to  the  WAC.  and  an 
area  of  sluggish  flow  between  this  southward  flow,  and  a 
flow  branching  from  the  Senj  gyre  northern  arm.  The 
northern  (upper)  ami  of  the  Senj  cyclonic  gyre,  observed 
and  model-generated  in  area  4  was  not  properly  resolved  in 
the  KO  model,  due  to  the  closeness  of  their  open  boundary. 
In  contrast  to  the  Trieste  gyre,  the  Senj  gyre  northern  ami  is 
unconstrained  geometrically,  enjoying  more  lateral  freedom 
and  exhibiting  less  rectilinear  (£  <  0.42)  flow. 

[47]  This  multigyral  response  is  a  prominent  generic, 
albeit  transient,  feature  of  the  northern  Adriatic  flow, 
imbedded  in  cyclonic  general  circulation.  As  several  previ¬ 
ous  modeling  (e.g.,  KO)  and  field  [e.g.,  Hopkins  et  al. , 
1999]  studies  have  established,  the  bora  generates  conver¬ 
gence  against  and  sea-level  rise  toward  the  western  bound¬ 
ary,  which  sustains  a  cyclonic  flow.  The  easterly  jets  of  both 
the  Trieste  and  Senj  bora  induce  flows  that  also  feed  the 
WAC.  Due  to  orographic  incisions  in  the  Dinaric  Alps,  the 
bora  exhibits  strong  horizontal  vorticity  (wind  stress  curl), 
which  is  responsible  for  imposing  the  observed  gyral 
pattern,  entailing  high-resolution  requirements  on  the  wind 
forcing  field.  Our  literature  review  reinforces  the  notion  that 
the  modeling  studies  based  on  coarser  ECMWF  wind  fields 
fail  to  capture  the  smaller-scale  vorticity  and  related 
response  pattern.  But  even  when  the  basic  response  pattern 
is  captured  correctly,  still  finer  resolution  is  likely  to  be 
required  to  resolve  the  important  details.  As  the  sea- 
response  circulation  pattern  is  tightly  coupled  w'ith  the 
forcing  wind  field,  any  weakness  in  wind  simulation  is 
promptly  reflected  in  the  modeled  sea  response. 

[4x]  The  COAMPS  model  successfully  simulated  the 
onset,  duration,  and  decay  of  the  wind  peaks  in  studied 
bora  episodes,  exhibiting  a  tendency  to  overpredict  the 
strength  of  the  bora  in  the  northern  Adriatic.  Bora  structures 
that  are  apparently  driven  by  topography  on  a  scale  finer 
than  resolved  by  COAMPS  (4  km)  are  not  well  represented 
in  the  simulated  winds,  so  improved  orography  representa¬ 
tion  may  improve  its  wind  field  predictions.  Further  im¬ 
provement  may  come  from  coupling  to  a  wave  model,  since 
feedback  of  wave-induced  stress  might  reduce  or  eliminate 
the  wind  speed  overprediction  observed  in  analyzed  bora 
episodes.  Using  a  coupled  atmosphere-ocean  wave  model, 
Doyle  [2002]  found  reduction  of  predicted  wind  speeds  in 
excess  of  10  ms  1  when  simulating  tropical  cyclones  and  a 
bora  event. 

[av]  The  bora  events  are  also  characterized  by  pro¬ 
nounced  cooling  and  heat  loss  to  the  atmosphere,  but  neither 
JRP  nor  WISE  had  a  component  to  adequately  address  this 
aspect  of  the  bora.  Our  NCOM  and  QUODDY  simulations 
were  therefore  focused  on  cooling  as  a  source  of  cold  water 


usable  as  a  tracer  of  bora-induced  flow  features.  The  model¬ 
ing  results  have  identified  the  shallow  NW  coastal  strip  as  an 
important  source  of  the  colder  water  observ  ed  in  a  sequence  of 
remotely  sensed  SST  fields  derived  from  AVHRR  data.  A 
comparison  with  a  related  sequence  of  NCOM-generatcd 
SST  fields  provided  good  semi-quantitative  agreement. 

[50]  With  ADCP-registered  current  time-series  in  excess 
of  200  days,  about  10  wave  and  tide-gauge  records.  5 
coastal  wind  data  time  series,  and  at  least  as  many  related 
series  generated  by  the  atmospheric  and  two  oceanic  mod¬ 
els,  there  is  great  potential  for  frequency  domain  analysis. 
Cross  and  cospectral  analyses  of  those  data  sets  should 
reveal  valuable  information  about  coherence,  energetics, 
and  phase  relationships  among  the  involved  variables, 
shedding  further  light  on  the  intricacies  of  the  bora-driven 
flows.  This  task  is  already  started  and  we  plan  to  submit  a 
comprehensive  report  in  a  separate  paper. 
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Science.  Education,  and  Sport  via  research  grant  0098 1 1 3  fhe  W  ISE 
Experiment  was  funded  by  the  Jadran  Project.  Insightful  comments  of  two 
anonymous  reviewers  helped  improve  the  manuscript.  The  AVHRR  images 
were  created  w  ith  a  graphic  interface  designed  and  programmed  by  Igor 
Tomaiic.  A  modified  vecDiagram  Matlab  function  from  Dave  Jones 
Fathom  Toolbox  was  used  to  create  the  PVDs.  COAMPS  is  a  registered 
trademark  of  the  Naval  Research  Laboratory  . 

References 

Barron.  C.  N..  A.  B  Kara.  P  J  Martin.  R.  C  Rhodes,  and  L.  I  Smedstad 
(2005).  Formulation,  implementation  and  examination  of  vertical  coordi¬ 
nate  choices  in  the  global  Navv  Coastal  Ocean  Model  (NCOM).  Ocean 
Model!.,  doi:  1 0. 1 0 1 6/j.ocemod.2005.0 1 .004. 

Beg  Paklar,  G.,  V.  Isakov.  D.  Koracin.  V  Kourafalou.  and  M  Orlic  (2001 ). 
A  case  studv  of  bora-driven  flow  and  density  changes  on  the  Adriatic 
Shelf  (January  1087).  Cont.  Shelf  Res.,  21,  1752  1783. 

Bergamasco.  A.,  and  M.  Gacic  ( 1996).  Barocltnic  response  of  the  Adriatic 
Sea  to  an  episode  of  Bora  wind. ./.  Phys.  Oi'eanogr.,  26.  1354  1369. 
Blumberg,  A.  F..  and  G.  L.  Mellor  (1987).  A  description  of  a  three- 
dimensional  coastal  ocean  circulation  model,  in  Three-Dimensional 
Coastal  Ocean  Models,  edited  by  N.  Heaps.  AGU.  Washington.  D.  C. 
Bogden,  P  S..  P  Malanotte-Rizzoli.  and  R.  Signell  (1996).  Open-ocean 
boundary  conditions  from  interior  data:  Local  and  remote  forcing  of 
Massachusetts  Bay.  J  Geophvs.  Res..  101.  6487  6500. 

Csanady.  G.  T.  (1975).  Hydrodynamics  of  the  large  lakes.  Ann  Rev  Fluid 
Mecli..  7,  357  386. 

Cushman-Roisin,  B..  M.  Gacic.  P.-M.  Poulain,  and  A.  Artegiani  (I  ds.) 
(2001).  Physical  Oceanography  of  the  Adriatic  Sea  Past.  Present  and 
Future.  304  pp..  Springer,  New  York 

Dorman.  C.  E..  et  al.  (2006).  February  2003  marine  atmospheric  conditions 
and  the  bora  over  the  northern  Adriatic  Sea.  J.  Geophvs  Res  . 
doi:  1 0. 1 029  2005 JC003 1 34,  in  press. 

Doyle.  J.  D  (2002).  Coupled  atmosphere-ocean  wave  simulations  under 
high  wind  conditions.  Mon.  Heather  Rev.  I.O),  3087  3099. 

Egbert.  G.  D..  and  S.  Y  Erofeeva  (2002).  Efficient  inverse  modeling  of 
barotropic  ocean  tides,  ./.  Atmos.  Oceanic  Techno!.,  IV.  183  204. 
Emery,  W.  J.,  and  R  E.  Thompson  (1997),  Data  Analysis  Methods  in 
Physical  Oceanography.  634  pp..  Elsevier.  New  York. 

Franco.  P.  and  A.  Michelato  (1992).  Northern  Adriatic  Sea  Oceanography 
of  the  basin  proper  and  of  the  w  estern  coastal  zone.  Si  1.  Total  F.nviron  . 
suppl..  35  62. 

Galperin.  B..  L  H.  Kantlia.  S  Hassid,  and  A  Rosati  (1988).  A  quasi¬ 
equilibrium  turbulent  energy  model  for  geophysical  flows.  ./  Atmos 
Sci..  45,  55  62. 


26  of  27 


C03SI3 


KUZMIC  ET  AL,:  NORTHERN  ADRIATIC  RESPONSE  TO  BORA  WIND 


C03SI3 


Hodur,  R,  M  (1997),  The  Naval  Research  Laboratory's  Coupled  Ocean'' 
Atmosphere  Mesoscale  Prediction  System  (COAMPS),  Mon,  Weather 
Rev..  125,  14(4  1430. 

Holland,  W.  Rr,  J.  C,  Chow,  and  F.  O.  Bryan  ( 1998).  Application  of  a  third- 
order  upwind  scheme  in  the  NCAR  Ocean  Model,  1  dim.,  11,  1487- 
1493. 

Hopkins,  T  S.,  A.  Artegiam,  C,  Kinder,  and  R  Parianie  (1999),  A  discus¬ 
sion  of  the  northern  Adriatic  circulation  and  flushing  as  determined  from 
the  ELNA  hydrography,  in  The  Adriatic  Sea,  Ecosyst.  Res.  Rep .  32, 
pp.  85  106,  Off,  for  Off.  Pub I.  of  the  Eur  Common,,  Luxembourg. 

Ivancan-Pieek,  B.,  and  V.  TuljS(|996),  A  case  study  of  severe  Adriatic  bora 
on  28  December  1992,  Tell  as.  Set:  .4,  48,  357  367. 

Janekovic.  I..  and  M.  Kuzmic  (2005),  Numerical  simulation  of  the  Adriatic 
Sea  principal  tidal  constituents,  Ann,  Geophvs.,  25,  3207  3218. 

Janekovic,  I.,  M.  Kuzmic.  J,  W.  Book,  and  H.  T.  Perkins  (2004).  A  tidal 
model  of  the  Adriatic  Sea:  ACT:  WISE  contribution  to  its  eurreni 
response  validation,  paper  presented  at  Venice  Adriatic  Workshop,  Off. 
of  Nav.  Res.  Global,  Venice,  Italy.  14  16  June. 

Kara.  A.  R  ...  C.  N.  Barron,  P.  J,  Martin.  L.  F.  Smedslad,  and  R.  C.  Rhodes 
(2005),  Validation  of  interannual  simulations  from  the  T8  Global  Navy 
Coastal  Ocean  Model  (NCOM),  Ocean  Model!.,  doi: I (1. 1016  j.ocemod. 
2005,01.003, 

Kondo.  J,  11975),  Air-sea  bulk  transfer  coefficients  in  diabatic  conditions. 
Boundary  Layer  MeteoroL,  9t  91  1 12, 

Kundu,  P.  K.  (1976),  Ekinan  veering  observed  near  the  ocean  bottom. 
V  Phys.  Oceanogr.,  6,  238  242. 

Kuzmic,  M.  (1991 ),  Exploring  the  effects  of  bora  over  the  northern  Adria¬ 
tic:  CZCS  imagery  and  a  mathematical  model  prediction,  hit.  J.  Remote 
Sens.,  12,  207  214. 

Kuzmic,  M..  and  M.  Orlic  ( 1987).  Wind-induced  vertical  shearing:  ALPEX 
M E  DA  LPEX  data  and  i  n ode!  I  i  ng  e  sere  i  se,  ,4  nn .  Geophvs. ,  58,  1 03  M2. 

Kuzmic.  M„  M.  Orlic,  and  Z,  Pasaric  {1988}.  Mathematical  modelling  (in 
Croatian),  in  1987  Annual  Report,  edited  by  R.  Prccali  and  Z.  Konrad, 
chap.  5,  pp.  331  352,  Joint  Yugoslav- Italian  Sci.  Coop.  Program, 
Rovinj,  Croatia 

Large.  W.  G.,  and  S.  Pond  (19811,  Open  ocean  momentum  flux  measure¬ 
ments  in  moderate  to  Strong  winds, ./.  Phys,  Oceanogr,,  II,  324  336, 

Lazie,  L..  and  I  ToSic  (1998),  A  real  data  simulation  of  the  Adriatic  Bora 
and  the  impact  of  mountain  height  on  the  Bora  trajectories.  Meteoml. 
Atmos.  Phys,,  66,  143  155. 

Lee,  C.,  er  al.  (2005a),  Nonhem  Adriatic  response  to  a  wintertime  bora 
wind  event,  Eos  Trans,  AGU .  <V6(16).  157  168, 

Lee,  C„  et  al.  (2005b),  Cruise  report:  DOLCE  VITA  I  and  2.  31  January 
23  February  and  26  May  15  June,  2003,  Tech.  Rep.,  APL-UW  TR  (1409. 

Loglisci.  N.,  M.  W.  Qian.  N,  Rachev,  C.  Cassardo,  A.  Longhetto.  R.  Turin i, 
P.  Trivero,  $.  Ferrarese.  and  C.  Giraud  (2004 1.  Development  of  an  atmo¬ 
sphere-ocean  coupled  model  and  its  application  over  the  Adriatic  Sea 
during  a  severe  weather  event  of  Bora  wind,  J.  Geophvs .  Res. ,  109, 
D0 1 102,  doi:  J0J  029/2003 JD003 95 6. 

Lynch,  D  R.,  and  F,  E.  Werner  (1991),  Three-dimensional  hydrodynamics 
on  Finite  dements.  Part  II:  Non-linear  time-stepping  model,  hit.  J.  Numer. 
Methods  Fluids ,  12,  507  533. 

Lynch.  D.  R.,  T.  C.  L  Justin.  C.  IT  Nannie.  and  F.  E.  Werner  (19961. 
Comprehensive  coastal  circulation  model  with  application  to  the  Gulf 
of  Maine.  Cant  Shelf  Res.,  16.  875  906. 

Martin.  P,  J.  (20001,  A  description  of  the  Naw  Coastal  Ocean  Model 
Version  L0,  NR1.  Rep.  NRUFR/7322-MU99621  42  pp.,  Nav.  Res.  Lab,, 
Sfennis  Space  Cent.,  Miss, 

Mel  tor,  G.  L.  ( 1991},  An  equation  of  slate  for  numerical  models  of  oceans 
and  estuaries. , l  Atmos.  Oceanic  Techno!.,  8,  609  611 

Mellon  G.  L.,  and  T.  Vamada  (1974),  A  hierarchy  of  turbulence  closure 
models  for  planetary  boundary  layers,  J .  Atmos.  Sci. ,  51,  1791  1806. 

Mellor,  G.  L„  and  T.  Vamada  1 1982),  Development  of  a  turbulence  closure 
model  for  geophysical  fluid  problems..  Rev.  Geophys.,  20,  85 1  875, 


Morey.  S.  L„  P.  J.  Martin,  J.  I.  O'Brien,  A.  A.  Wall  craft,  and  J.  Zavala- 
Hidalgo  (2003),  Export  pathways  for  river  discharged  fresh  water  in  the 
northern  Gulf  of  Mexico,  J.  Geophvs.  Res..  1081  CIO).  3303,  doi:  10. 1029 
2OO2JC00 1 674. 

Orlic,  M.t  M.  Kuzmic,  and  Z.  Pasaric  (1994  k  Response  of  the  Adriatic  Sea 
to  the  bora  and  sirocco  forcing.  Cant ,  Shelf  Res .  14,  91  116. 

Penzar,  B.  (1977),  Air  pressure  wind  (in  Croatian),  Ptihti  [mzmivanju 
vremena  i  fell  me.  2.  117  pp. 

Perkins,  H.  T.,  F.  de  Strobe  I,  and  L.  Gualdesi  (2000),  The  Barney  Sentinel 
Trawl- resistant  ADC  P  bottom  mount:  Design,  testing,  and  application. 
IEEE  J.  Oceanic  Eng.,  25,  430  436. 

Poulain,  P.-M.,  V.  H.  Koumfalou,  and  B  Cushmun-Roisin  (2(101 ),  Northern 
Adriatic  Sea.  in  Physical  Oceanography  of  the  Adriatic  Sea:  Past.  Pre¬ 
sent  and  Future,  chap.  5.  edited  by  B.  Cushman-Roisin  et  al,T  pp  143 
165,  Springer.  New  York. 

Pullen,  J.,  J.  D.  Doyle,  R.  Hodur,  A.  Ogston.  J.  W  Book,  H  Perkins,  and 
R.  Signell  (2003),  Coupled  ocean-atmosphere  nested  modeling  of  the 
Adriatic  Sea  during  winter  and  spring  2001,  J.  Geophvs  Res., 
I081C  i  0).  3320.  doi:  1 0,  J  029  2003 JCO0 1 780. 

Rachev,  N.,  and  R  Purini  (20(11  ),  The  Adriatic  response  to  Ihe  bora  forcing 
A  numerical  study,  Nuovo  Cimento,  24 C.  303  31 1 . 

Raicich.  R.  (1994).  Note  on  the  flow  rates  of  the  Adriatic  rivers.  Tech  Rep. 
RF  02  94,  8  pp..  Cons.  Naz,  delta  Ric.  Inst.  Sper.  Talassogratlco,  Trieste, 
Italy. 

Schwab,  D  J.,  A  H,  elites.  C  R.  Murthy.  J,  E.  Sand  well,  L.  A.  Meadows, 
and  G  .  M.  Meadow  s  (1989),  The  effect  of  wind  on  transport  and  circula¬ 
tion  in  Lake  St,  Clair,  J.  Geophys.  Res,,  94,  4947  4958, 

Sc h wing,  F.  B.,  L.-Y.  Oey,  and  J.  O.  Blanton  (1985),  Frictional  response  of 
continental  shelf  waters  to  local  wind  forcing,  J .  Phvs.  Oceanogr ,  15, 
1733  1 746. 

Shtokman.  V.  R.  (1941).  Vetrovoi  nagon  i  gorizontalynaya  cirkulaciya  v 
zamknutom  morye  ncbolshoy  glubiny  (in  Russian),  hr.  ,4V  SSSR.  Set: 
Geografk  h es ka\  a  i  geofiricheskaya 

Signell.  R.  P.  S.  Camiel,  L.  Cavaleri,  J.  Chiggiato,  J.  D.  Doyle.  J,  Pullen, 
and  Mr  Sclav  o  (2005).  Assessment  of  wind  quality  for  oceanographic 
modelling  in  semi-enclosed  seas, ,/.  Mar.  Sysr.,  55,  217  233. 

Smagorinsky.  J.  ( 1963).  General  circulation  experiments  with  the  primitive 
equations  I.  The  basic  experiment,  Mon.  Weather  Rev,  91,  99  164. 

Spoler-Canic.  K.,  and  L  Kraljevic  (2005),  The  impact  of  model  resolution 
on  simulated  distribution  of  hura,  Croat.  Meteoml  J.,  40.  296  299, 

Stravisi,  F.  (1977),  Bora  driven  circulation  in  the  northern  Adriatic,  Bolt. 
Geode  Tear,  Appl. ,  19,  95  102, 

Thompson,  K.  R  and  D.  I.  Pugh  1 1986).  The  subudal  behaviour  of  the 
Celtic  Sea  II.  Currents,  Co  tit.  Shelf  Res..  5 ,  321  346. 

Wang,  X,  H,  (2005),  Circulation  and  heat  budget  of  the  northern  Adriatic 
Sea  (Italy)  due  to  a  Bora  event  in  January  2001:  A  numerical  model 
study.  Ocean  Model!.,  10,253  271. 

Yoshino,  M,  ( 1972).  Monthly  mean  wind  frequency  and  force  by  directions 
ar  selected  climatological  stations  in  Yugoslavia.  Climatai.  Votes,  10, 
23  44. 

Yoshino,  M,  (1976),  Local  Wind  Bora,  289  pp.,  Univ.  of  Tokyo  Press, 
Tokyo, Japan. 

Zore-Armanda,  M.,  and  M,  Gacie  (1987),  Effects  of  bora  on  the  circulation 
in  the  North  Adriatic,  Ami.  Geophvs,,  58,  93  102, 


J.  Wr  Rook  and  P.  J,  Martin.  Oceanography  Division.  Naval  Research 
Laboratory,  S tennis  Space  Center,  MS  39529^1 1 05,  USA. 

J  D,  Do  vie.  Marine  Meteorology  Division,  Naval  Research  Laboratory. 
Monterey,  CA  93943.  USA. 

I.  Janekovic  and  M,  Kuzmic,  Center  tor  Maune  and  Environmental 
Research,  Ruder  BoSkovic  institute.  P.O.  Box  180,  100(12  Zagreb.  C  roatia. 
( kuzm  icfrr  nidjer.  i  rb.  hr ) 


27  of  27 


